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FOREWORD 


This  work  has  been  conducted  to  expand  upon  previous  work  and  to  provide  a  design  tool  for 
use  in  estimating  the  aerodynamics  of  today’s  high  performance  tactical  weapons.  The  resulting 
computer  program,  into  which  this  work  has  been  incorporated,  allows  one  to  predict  performance 
and  to  conduct  a  static  and  dynamic  stability  analysis  in  the  preliminary  and  intermediate  design 
stages  without  costly  and  time-consuming  wind  tunnel  tests. 

Support  for  the  work  was  provided  by  the  following  sponsors: 


The  Naval  Sea  Systems  Command  under  the  Surface-Launched  Weapons  Aerodynamics  and 
Structures  Block  Program 
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Technology  Block 

The  U.S.  Army  Missile  Command  under  Project  No.  1L162303A214 
The  Office  of  Naval  Research  under  Project  No.  F4141 1 


Major  procurement  contracts  from  the  above  funds  were  let  to  North  Carolina  State  Univer¬ 
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This  report  presents  the  theoretical  basis  for  the  fourth  increment  in  the  development  of  a 
computer  code  for  the  rapid  prediction  of  aerodynamics  of  tactical  weapons  in  the  region 
0<M^<8and0<a<  180°. 

The  general  approach  has  been  to  combine  existing  and  newly  developed  approximate  com¬ 
putational  methods  into  a  single  computer  program  to  compute  aerodynamics.  Computational 
times,  required  for  the  estimate  of  static  and  dynamic  aerodynamic  coefficients  for  a  body-tail- 
canard  configuration  for  one  freestream  condition,  are  in  CPU  seconds  on  the  CDC  6700  as  opposed 
to  minutes  or  hours  required  for  more  detailed  physical  and  numerical  models.  The  accuracy 
obtained,  however,  is  compatible  with  that  required  for  preliminary  or  intermediate  design 
estimates. 

The  first  increment*  led  to  the  development  of  a  code  for  the  prediction  of  static  aerodynamic 
coefficients  to  Mach  number  3  and  angles  of  attack  to  15°  for  body-alone  configurations. 

A  second  increment^  supplemented  the  earlier  work  and  resulted  in  a  code  for  the  prediction 
of  static  aerodynamic  coefficients  for  body  alone,  body-tail,  and  body-tail-canard  configurations  for 
the  same  range  of  freestream  conditions. 

A  third  increment^  supplemented  the  work  of  the  first  two  increments  and  extended  the 
computational  capability  to  include  the  prediction  of  dynamic  derivatives,  again  for  the  same  range 
of  freestream  conditions. 

The  prime  objective  of  the  current  work  effort  is  to  extend  the  capability  of  the  aeroprediction 
code  to  the  higher  Mach  number  and  higher  angle-of-attack  region  of  advanced  designs  (M^^  -*■  8, 
a  -*■  180°).  A  secondary  objective  is  to  improve  the  estimate  of  static  derivatives  in  the  transonic 
Mach  number  regime  and  to  improve  the  dynamic  derivative  prediction  capability  at  all  Mach 
numbers. 


ANALYSIS 


CONFIGURATIONAL  GEOMETRY 

The  basic  configurational  geometry  remains  the  same  as  for  References  2  and  3.  For  various 
Mach  number  and  angle-of-attack  regions  there  are  geometric  restrictions.  These  will  be  elaborated 
upon  in  later  sections. 

Intended  applications  are  for  spin-stabilized  projectiles,  bombs,  and  rockets  (after  burn-out)  in 
the  unguided  class  and  guided  missiles/projectiles.  Inlet  and  plume  effects  for  various  weapons  are 
not  considered.  The  code,  however,  has  on  occasion  been  used  for  some  nontactical  configurations. 


The  degree  of  the  configurational  complexity  that  can  be  considered  is  illustrated  in  Figure  1. 
This  most  complex  body-of-revolution  consists  of  a  spherical  nose  cap.  two  piecewise-continuous 
nose  sections,  a  straight  afterbody,  and  a  boattail.  The  wing  or  canard  has  a  trapezoidal  planform 
with  a  biconvex  or  modified  double  wedge  cross  section  and  sharp  or  spherically  blunted  leading  or 
trailing  edges.  Tip  edges  are  assumed  parallel  to  the  freestream.  Twist  or  airfoil  distortion  is 
neglected  and  piecewise  similar  airfoil  shape  variation  with  span  is  assumed. 

REVIEW  OF  PREVIOUS  METHODS 
Body-Alone  Static  Aerodynamics 

A  summary  of  the  various  methods  for  computing  body-alone  aerodynamics  appears  in 
Figure  2.  The  majority  of  the  methods  used  were  adapted  from  those  available  in  the  standard 
literature.'*'^  The  empirical  and  semiempirical  schemes  used  for  the  transonic  lift  and  wave  drag 
are  presented  in  some  detail  in  Reference  1.  The  combined  Newtonian-perturbation  theory,  also 
presented  in  Reference  1 ,  was  developed  so  that  reasonable  results  for  static  aerodynamics  could 
be  obtained  at  low  supersonic  Mach  numbers  for  blunt-nosed  configurations.  Mach  numbers  0.8 
and  1.2  are  normal  division  points  for  the  three  Mach  number  regions. 

Wing  and  Interference  Static  Aerodynamics 

The  methods  used  to  compute  lifting  surface  alone  and  interference  aerodynamics  are  shown 
in  Figure  3.  Methods  adapted  from  the  standard  literature  were  taken  from  References  10  through 
16.  The  remaining  methods  for  wing-body  interference,  trailing  edge  separation  drag,  body  base 
pressure  caused  by  tail  fins,  and  wing  drag  numerical  techniques  are  detailed  in  Reference  2. 


Figure  1 .  General  Configurational  Geometry 
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Dynamic  Derivative  Computational  Methods 


Methods  for  computing  the  dynamic  derivatives  are  listed  in  Figure  4.  Applicable  references 
taken  from  Reference  3  are  10  through  12  and  15  through  26. 

HIGH  MACH  NUMBER  STATIC  AERODYNAMICS 

This  work  was  performed  under  contract  by  F.  De  Jarnette  of  North  Carolina  State  University 
during  FY  77-FY  80. 

No  real  gas  effects  are  considered,  although  some  real  gas  effects  become  important  above 
M^  =  5.  Strong  viscid-inviscid  interactions  are  neglected.  Lifting  surface-body  interference  is 
neglected.  The  body  base  pressure,  wing  base  pressure,  body  and  wing  friction  routines  are  the 
same  as  those  reported  in  References  1  and  2. 

One  consideration  is  the  Mach  number  range  of  validity  of  the  high  Mach  number  prediction 
method.  This  will  be  discussed  in  more  detail  in  the  methods  evaluation  section.  A  Mach  number, 
M  =  M^j,  will  be  the  dividing  Mach  number  at  the  lower  range  limit  and  is  a  code  input. 

Body-Alone  Inviscid  Aerodynamics 


The  method  used  for  pointed  bodies  with  attached  shock  waves  is  based  on  a  modification  of 
the  method  given  in  Reference  27.  For  blunt  bodies,  a  modified  Newtonian  pressure  distribution  is 
matched  to  the  second-order  shock-expansion  pressure  distribution.  This  method  is  a  modification 
of  that  used  in  Reference  28. 
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Figure  4.  Methods  Used  to  Compute  Dynamic  Derivatives 
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The  perturbation  theory  of  Reference  4  breaks  down  when  the  local  slope  is  greater  than  the 
Mach  angle.  For  more  slender  bodies,  other  nonlinear  terms,  which  are  neglected  by  perturbation 
theory,  become  significant  at  higher  Mach  numbers  even  though  the  local  slopes  are  smaller  than 
the  Mach  angle.  Modified  linear  theory  of  the  type  presented  in  Reference  29  uses  the  local  P  in  an 
iterative  solution  and  avoids  the  slopc-Mach  angle  limitation.  The  approach  of  Reference  4  is 
probably  not  applicable  at  Mach  numbers  above  4.  The  current  approach  is  given  in  Reference  30. 


For  a  body  at  small  angles  of  attack,  the  deviation  of  free  streamlines  from  meridian  lines  is 
small  and  hence  is  neglected. 


From  the  method  of  characteristics,  for 


^  _  27P  ^ 

3s  sin  2  /./  3s 

ap 

ac, 


a  body  of  revolution  on  a  free  streamline 

_! _ ^ 

cos /a  3C| 

27P  /  35  sin  p  sin  S\ 

sin  2  tx  Uc,  ^  r  / 


(1) 

(2) 


Cl  refers  to  the  left  running  or  outgoing  characteristic.  Incoming  characteristics  are  neglected 
(i.e.,  reflection  from  a  bow  shock).  The  angle,  5,  is  the  angle  of  the  local  streamline  tangent;  s  is 
the  distance  along  a  streamline.  The  original  body  is  replaced  by  a  tangent  body  as  shown  in 
Figure  5. 


Figure  5.  Tangent  Body  Geometry 


The  solution  starts  as  a  cone  solution  or  as  a  solution  that  matches  the  Newtonian  pressure 
distribution.  The  match-point  procedure  will  be  discussed  in  detail  in  a  later  section. 

At  a  comer  such  as  1-2,  the  terms  on  the  right-hand  side  of  Equation  1  are  neglected  com¬ 
pared  with  the  left-hand  side.  Thus,  at  the  comer,  the  well  known  Prandtl-Meyer  relationship 
holds. 

On  the  straight  section  2-3,  35/ds  =  0.  If  one  considers  a  free  streamline  close  to  the 
surface  2-3,  and  assumes  that  the  pressure  is  constant  on  that  streamline,  that  the  Cj  characteristics 
are  straight,  that  the  surface  Mach  number  is  a  constant,  and  that  AC]  is  a  constant,  then  one  can 
derive  approximately  (at  a  =  0  meridian  and  streamlines  coincide) 

p  =  k,  +  kje^J®'.  (3) 


Here,  s'  =  x'/cos  5  is  the  distance  from  the  corner  1-2,  x'is  the  corresponding  horizontal  distance 
from  the  comer,  and  kj  =  p^  is  the  asymptote  for  the  pressure  distribution.  For  6  >  0  on  Section  2-3 
and  s'  very  long,  p^  equals  the  cone  pressure  at  the  freestream  Mach  number.  For  6  <  0,  the  assump¬ 
tion  is  made  that  p^  =  p^o  or  the  freestream  pressure.  This  last  assumption  can  be  relatively  poor. 
For  5  <  0,  a  pressure  distribution  approximating  the  recompression  on  an  conical  boattail  is 
is  needed.  In  general.  Equation  3  must  be  considered  heuristic.  The  remaining  constants  can  be 
determined  from  a  knowledge  of  the  pressure,  P2,  and  the  pressure  gradient,  (dp/dslj. 

The  pressure  distribution  is  then  given  by 


P  =  Pa  +  (P2  -  Pa)®”’' 


\  va/2 

(P2  -  Pa) 


(4) 

(5) 


where  P2  is  obtained  from  the  Prandtl-Meyer  relationship  for  an  isentropic  expansion.  Also, 

^2  -  1^1  =  -(82  -  5] )  (6) 

and 


=  -/y  =  /p-]  /fri  -  tan->  v/m2  -  1 


(7) 
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where 


X  = 


U  = 


2yp 

sin  2ju 

-1  1 
sin  >  ^ 


7 


(8) 

(9) 


(10) 


and  Po  is  the  stagnation  pressure  behind  a  bow  shock. 

For  a  compression  comer,  the  isentropic  relationship  is  used.  If  an  isentropic  solution  for 
P2  is  not  obtainable,  p  =  p^  is  assumed  for  the  straight  section.  The  term  Op/ds)2  is  obtained  from 
the  approximate  isentropic  relationship  from  Reference  27.  Thus 


where 


and 


7pM^ 

”  2(M2  -  1) 


7+1 

T~ 


7-H 

2(7-1) 


(12) 


(13) 


The  same  relationship  is  used  for  compression  comers. 

The  pressure  relationship  breaks  down  when  rj  <  0  and  an  asymptote  cannot  be  reached. 
When  Tj  <  0,  p  =  Pj  for  an  expansion  comer  and  p  =  P2  for  a  compression  comer  are  chosen. 

The  pressure,  which  is  used  for  load  integration,  is  evaluated  at  the  tangency  points  (C)  in 
Figure  5. 

Equation  4  provides  a  solution  along  a  meridian  plane  surface  streamline.  For  small  angles  of 
attack,  a  pressure  coefficient  expansion  is  assumed  of  the  form 


(14) 


Cp  =  (Cp)t,^Q  -  A  sin  2a  cos  0  +  (F  cos^i^  +  A  sin^^)  sin^a/Z 
=  Cp,  +  (Cp2  -Cpje-" 

where  (p  is  the  azimuth  angle  measured  from  the  leeward  plane 


(IS) 

(16) 

and 

= C:) 

\  3a2  4^o,0  =  jr/2 

(17) 

The  streamlines  make  an  angle,  e,  with  the  meridian  lines  of  the  order  e  ~  sin  0  sin  a.  s'  is 
along  a  surface  streamline.  Furthermore,  P2  and  its  a  derivatives  are  independent  of  0  since  e  is 
small;  p,,  however  is  dependent  on  0. 

It  may  be  readily  shown  that  the  loading  functions  (Cp)a=o>  C,  and  A  have  the  same 

functional  form  as  Equation  4;  i.e. 

(Cp)a  =  0  “  (Cpa)a  =  0  (Cp2  ~  Cpa  )(J;  =  0  ^  ^ 

(18) 

A  =  A,  +  (A2  -  Aj)  e~’^ 

(19) 

r  =  r,  +  (r2-r,)  e-'i 

(20) 

A  =  A,  +  (A2  -Aa)  e-’J 

(21) 

Cp2  is  given  in  terms  of  Cpi  from  the  Prandtl-Meyer  relationships.  Differentiation  of  the 
Prandtl-Meyer  relations  with  respect  to  a  once  and  then  twice  leads  to  the  following  relationships: 


A2 


(22) 


^2 


(23) 


^2 


h 


+ 


27MiAf 


[G(M2)-G(M,)1  \ 


(24) 
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where 


and 


G(M)  = 


M**  -  2M^  +  2 


i  -  n3/2 


(M-'  -  1) 


(25) 


(26) 


The  expressions  for  (Cpj)^  =  0’  ^a>  f'a’  obtained  from  an  approximation  for  the 

asymptotic  pressure  coefficient;  i.e.. 


Cpa  ^^pa^O  =  0  ^^p 


(27) 


where  (Cpj,)a  =  o  2®’’°  for  fi  <  0  and  for  6  >  0  is  obtained  from  an  approximate  relation  for  the  cone 
solution  taken  from  Reference  3 1  where 


where 


(Cpa)a  =  ( 


sin^6 


|,  +  2 
[  (7-  DK^  +  2 


8n 


7  +  > 


=  (M^  -  Usings 


(28) 


(29) 


Equation  28  is  accurate  except  upon  approaching  the  detachment  Mach  number. 


Next,  ACp,  the  angle-of-attack  dependence  of  Cp^.  is  given  by  a  blend  of  slender  body  and 
Newtonian  approximations;  i.e. 


AC_  =  -sin  2a  sin  26  cos0  +  sin^acos2  6 


(1  -  tan^  6)  -  (2  + 


t) 


(30) 


From  Equations  15  through  17  and  30  one  can  obtain 
Aj  =  sin  25 

r,  =  2  cos^6 


tan^6) 


(31) 

(32) 
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and 


A,  =  :  cos-5 


(=-|) -(2+1) 


(33) 


In  order  to  start  the  solution,  Cp, ,  Aj ,  A, ,  and  P,  must  be  specified  as  initial  values.  For  a 
pointed  or  truncated  body,  these  relations  are  given  by  Equations  28.  29,  and  31  through  33.  Fora 
blunt  body  (not  necessarily  spherical),  the  modified  Newtonian  pressure  distribution  may  be 
matched  to  Equation  14  at  some  match  point.  The  match  point  is  taken  as  M  =  1.1,  as  determined 
from  the  isentropic  pressure  relationship  combined  with  the  modified  Newtonian  pressure  relation¬ 
ship  evaluated  at  a  =  0. 


The  modified  Newtonian  pressure  distribution  is  given  as 

~  ^po  (‘■"O!' «  sin  6  -  sin  a  cos  6  cos  0)^  (34) 

Comparing  Equations  14  and  34  yields 

<^pi  =  (^po  (35) 

A,  =  Cpu  sin  5  cos 6  (36) 

P,  =  2Cp^  cos  2  6  (37) 

and 

A,  =  -2Cp„  sin^S  (38) 

Evaluation  of  the  loading  function  Equations  18  through  21  at  the  end  of  an  interval  (point 
3  in  Figure  5)  continues  the  solution. 

Integration  of  the  pressure  distribution  yields  the  loading  aerodynamic  coefficients  as 


f  RR’  Wp)a-o 

.  P  +  A  .  2  1  . 

+  — ^ —  sin^a  dx 

- 

(39) 


and 


4  sin  2a 


RA  dx 


C 


m(about  nose) 


-4  sin  2a 


I  RA(x  +  RR')dx 
'b 


(40) 


(41) 


All  dimensions  above  are  in  calibers.  The  reference  area  coincides  with  the  maximum  nose 
radius.  The  second  term  in  Equation  41  is  neglected. 
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Wing-Alone  and  Interference  Aerodynamics 


The  geometry  used  for  the  wing  is  the  same  as  that  used  by  Moore.^  Two  basic  airfoil  shapes 
are  used-the  double  wedge  and  the  biconvex.  Each  of  these  airfoil  shapes  may  have  a  sharp  or 
blunt  leading  or  trailing  edge.  Figure  6  shows  the  planform  and  airfoil  sections  for  the  double 
wedge,  and  Figure  7  shows  them  for  the  biconvex  case.  No  twist  or  airfoil  distortion  is  allowed. 
Piecewise  similar  airfoil  shape  variation  with  span  is  assumed. 

The  theory  used  to  determine  the  pressure  distribution  is  based  on  two-dimensional  supersonic 
flow  properties.  At  a  given  point  on  the  surface,  the  slope  of  the  surface  relative  to  the  undisturbed 
freestream  is  determined  first.  If  this  angle  makes  a  compression  surface,  the  pressure  is  calculated 
from  tangent  wedge  theory.  In  this  theory,  the  angle  between  the  surface  and  the  freestream  is 
the  wedge  angle  used  to  calculate  the  pressure.  If  the  angle  indicates  an  expansion  surface,  the 
pressure  is  calculated  from  Prandtl-Meyer  theory  for  an  expansion  from  freestream  properties 
through  this  angle.  These  angles  are  illustrated  in  Figure  8,  and  the  two  theories  are  referred  to  as 
the  tangent-wedge  and  Prandtl-Meyer  theories.  For  blunted  leading  edges,  modified  Newtonian 
theory  is  used. 
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SECTION  F  F  (SHARP  LEADING  AND  TRAILING  EDGES) 
Z 


SECTION  F  F  ISLUNT  LEADING  AND  TRAILING  EDGES) 
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Figure  6.  Wing  with  Modified  Double  Wedge  Airfoil  Section 


SECTION  F  F  (BLUNT  LEADING  AND  TRAILING  EDGES! 
Z 


Figure  7.  Wing  with  Biconvex  Airfoil  Section 
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PRANDTL  MEYER 
EXPANSION  <  01 


Figure  8.  Wing  Pressure  Distribution 


Here 

fi  =  unit  vector  normal  to  the  surface 

A 

=  unit  freestream  vector  and 

6eff  =  effective  deflection  angle  =  sin"i  (-n  •  V^) 

For  a  wing  connected  to  a  body,  the  two-dimensional  theories  described  above  are  in  error 
near  the  body  and  the  tip  of  the  wing.  However,  Hilton^^  showed  that  for  a  rectangular  wing  these 
errors  cancel  each  other  (see  Figure  9).  For  wings  other  than  rectangular  ones,  the  errors  do  not 
exactly  cancel,  but  the  predicted  forces  and  moments  have  been  found  to  be  reasonably  accurate. 

If  the  wing  pressure  is  calculated  by  two-dimensional  theory  (shock-expansion): 

1.  The  pressure  near  the  body  is  too  low  because  of  wing-body  interference 

2.  The  pressure  in  the  tip  region  is  too  high  because  of  tip  losses 

These  effects  cancel  each  other  and  the  two-dimensional  theory  gives  the  correct  force  on  a 
rectangular  wing. 

The  body  ahead  of  the  most  forward  wing-body  junction  point  affects  the  wing  loading  as 
does  the  body  close  to  the  root  chord  junction  line.  The  wing  affects  the  body  loading  downstream 
from  the  forward  wing-body  junction  point.  Slender  wing-body  interference  approximations  used 
in  Reference  2  are  valid  only  for  the  low  Mach  number  flow  regime. 
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Figure  9.  Wing-Body  Interference 


The  double-wedge  wing  shown  in  Figure  6  has  three  flat  panels  on  the  top  and  three  flat  panels 
on  the  bottom.  The  pressure  on  each  panel  is  determined  by  the  tangent  wedge  or  Prandtl-Meyer 
theory,  and  the  pressure  is  constant  over  each  panel.  If  the  leading  edge  is  blunted,  modified 
Newtonian  theory  is  used  to  calculate  the  pressure  over  the  blunted  portion  in  the  same  manner 
used  by  Moore. ^  If  the  trailing  edge  is  blunted,  the  empirical  two-dimensional  base  pressures  given 
by  Moore^  are  used  and  then  extrapolated  for  the  higher  Mach  numbers. 


Each  flat  panel  is  represented  by  the  x,  y.  and  z  coordinates  of  the  four  comer  points  of  the 
trapezoidal  surface  as  shown  in  Figure  10. 

Referring  to  Figure  10,  V|_4  is  the  vector  connecting  points  1  to  4  and  V2-3  connects  2  to  3. 

Then 

V|_4  =  Ax4r  +  Ay4j'  +  Az4k  (42) 

A  A 

where^i,  j,  and  k  are  the  unit  vectors  in  the  x,  y,  and  z  directions,  respectively,  and 


and 


AX4  =  X4  -  X| 

Ay4  =  y4  -  Vi  =  b/2  (43) 

AZ4  ~  Z4  “  Z| 

V2_3  =  Ax3i  +  Ay3j  +  Az3k  (44) 
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Figure  10.  Geometry  of  Wing  Panel 


where 

Axj  =  X3  -  X2 

=  Ya  “  Ya  =  b/2  (45) 

AZ3  =23-22 

The  unit  vector  normal  (outer)  to  the  panel  is  then 

.  V,_4XV2-a 

n  =  — -  (46) 

IV,_4  X  V2-3I 

Relative  to  the  wing  coordinate  system,  the  unit  vector  parallel  to  the  freestream  velocity  vector  is 

V„  ~  cos(a  +  6)i  +  sin(a  +  6)j  (47) 

where  o  is  the  angle  of  attack  and  5  is  the  deflection  angle  of  the  wing.  The  effective  angle  between 
V„  and  n  is  called 

where 


sin  =  -n  • 

Performing  the  indicated  vector  operations  with  the  equations  above,  one  obtains 

[cos  (a  +  6)(A24  -  A23)  +  sin  (O' +  6)(Ax3  -  Ax4)lb/2 


sin  6jff  = 


D 


(48) 


(49) 
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where 


D  =  {[b/2(Az3  -  Az4)|-  +  [AX3AZ4  -Ax4Az3)^  +  [b/2(Ax4  -  Axj)]^}'''-  (50) 

If  5(.„  >  0,  the  tangent  wedge  theory  is  used  with  a  wedge  angle  of  5„,f.  The  explicit  wedge  solution 
given  by  Reference  34  is  used  here.  If  5^,,-  <  0,  the  Prandtl-Meyer  expansion  is  used  with  |6gff|  as 
the  expansion  angle  from  freestream  conditions. 

Once  the  pressure  on  each  panel  is  calculated,  the  axial  force,  normal  force,  and  pitching 
moment  can  be  calculated.  The  effective  area  for  the  axial  force  on  a  panel  is 

Sa  =  1^2  “  Zi  Z4  -z3]b/4  (51) 

and  therefore  the  axial  force  is 

A  =  pSa  (52) 

where  p  is  the  pressure  on  the  panel.  Note  that  Sa  may  be  positive  or  negative.  The  effective  area 
for  the  normal  force  on  a  panel  is 

Sn  =  [^2  “  Xi  +  X4  -  X3lb/4  (53) 

and  the  normal  force  is 

N  =  ±pSn  (54) 

where  the  +  sign  is  used  for  panels  on  the  lower  surface  and  the  -  sign  for  panels  on  the  upper 

surface  of  the  wing. 

For  biconvex  airfoils,  it  is  convenient  to  divide  the  wing  into  a  front  half  for  x  <  C/2  and  a 
rear  half  for  x  >  C/2.  The  chordwise  pressure  distribution  is  calculated  for  the  chord  at  midspan. 
The  geometry  of  the  wing  is  shown  in  Figure  1 1. 

The  airfoil  section  for  the  front  half  of  the  wing  is  considered  first  as  a  cross  section  perpendic¬ 
ular  to  the  leading  edge.  In  this  cross  section  the  airfoil  is  represented  by  a  circular  arc  from  the 
midchord  position  forward  to  the  nose.  If  the  leading  edge  is  sharp,  this  circular  arc  extends  to  the 
leading  edge  itself  However,  if  the  leading  edge  is  blunted,  the  circular  arc  blends  into  the  cylin¬ 
drical  leading  edge.  From  the  geometry  of  similar  triangles,  the  slope  of  the  circular  arc  at  its  most 
forward  position  is  given  by 

-2(cos  ALi  Cavg/^-fLi  )  (ru:  - 

sin  6,  =  - (55) 

(cosAL,.;Cavp/2-r,  |.;)2  +  (rL|.  - 
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AIRFOIL  OF  FRONT  HALF  AT  MIDSPAN 
(X’,  Y',  Z  COORDINATES) 

Figure  1 1 .  Geometry  for  a  Biconvex  Airfoil 


and  the  coordinates  of  this  point  are 


X|  =  rL|;(l -sin  6|)/cos  Alj; 

(56) 


Z,  =  Fli.;  cos  5, 

Now  consider  the  airfoil  section  in  a  plane  perpendicular  to  the  y-axis.  In  this  plane  the  circular 
arc  described  above  is  approximated  by  the  parabola 


z 


tan  5 1  cos  Alj; 


(57) 


This  expression  can  be  used  to  determine  the  effective  angle  for  the  tangent  wedge  or  Prandtl- 
Meyer  theories  as 


avg 


tan6(  cos  Ale  - cos  (a +  6) 


sin  5(,,c  = 


C'iivg 
X, - — 


±  sin  (a  +  6) 


1  + 


(^)j 


21 1/2 


(58) 


where  the  +  sign  corresponds  to  the  lower  surface  and  the  -  sign  to  the  upper  surface.  For  the 
front  half  of  the  wing 


Xl 


(59) 


The  pressure  is  calculated  at  10  positions  along  this  arc  and  integrated  by  Simpson’s  rule  to  obtain 
the  axial  force,  normal  force,  and  pitching  moment. 

The  rear  half  of  the  wing  is  handled  in  the  same  fashion  as  the  front  half. 

For  configurations  with  wings  and  canards,  the  analysis  described  above  is  used  to  treat  each 
configuration  like  a  separate  wing.  The  contributions  for  each  component  are  then  combined  to 
determine  the  overall  forces  and  moment. 

HIGH  ANGLE-OF-ATTACK  AERODYNAMICS 

Two  basic  approaches  have  been  used  to  predict  high  angle-of-attack  aerodynamics. 

The  first  of  these  is  an  empirical  approach. 25  Body-alone,  fin-alone,  fin-in-the-presence-of-a- 
body,  and  interference  aerodynamic  coefficient  components  are  obtained  from  empirical  functional 
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form  fits  for  a  family  of  body-tails  for  a  given  range  of  angles  of  attack,  Mach  numbers,  and  roll 
angles.  No  control  deflections  are  considered. 


The  second  approach  is  much  more  elaborate.^  It  uses  the  semiempirical  or  heuristic  cross- 
flow  analogy  with  wing-alone  loading  and  other  data  for  body  vortex  parameters.  Vortices  from 
the  nose,  body,  and  lifting  surfaces  are  tracked  in  a  stepwise  manner  with  distance  along  the  body 
centerline.  Control  deflections  are  considered.  Roll  moment  and  side  moment  coefficients  are 
predicted  as  well  as  longitudinal  aerodynamic  coefficients.  Body-canard-tail  configurations  as  well 
as  body-tail  configurations  are  considered. 

Currently,  the  aeroprediction  code  contains  the  coding  from  References  35  and  37.  The 
method  is  restricted  to  body-tail  configurations.  However,  it  requires  about  one-tenth  of  the 
storage  required  for  the  code  in  Reference  36,  and  computational  times  are  much  faster.  For  the 
currently  used  method,  the  range  of  input  parameters  for  a  cruciform  configuration  are 

1.  Mach  number:  0.8  to  3.0 

2.  Angle  of  attack;  0°  to  180°  for  isolated  components  (roll  angle  =  0°)  and  0°  to  45°  for 
body-tail  combinations  at  arbitrary  roll  angles  from  0°  to  1 80° 

3.  Tail:  Trapezoidal  planform,  edges  parallel  to  body  centerline 

a.  Leading  edge  sweep  angle:  0°  to  70° 

b.  Taper  ratio:  0  to  1 

c.  Aspect  ratio  (two  fins);  0.5  to  2.0 

4.  Nose  length  (pointed  tangent  ogive);  1.5  to  3. „  calibers 

5.  Cylindrical  afterbody;  6  to  18  calibers  length  and  with  the  end  of  the  body  parallel  to 
the  tail  trailing  edge 

6.  Total  span-to-diameter  ratio  (two  fins):  1  to  3-1/3. 

Roll  angle  definition  and  positive  fin  load  orientation  is  shown  in  Figure  12.  Note  that  0  is 
measured  180°  from  the  definition  in  the  aeroprediction  code  (i.e.,  the  windward  side  instead  of 
the  leeward  side). 

The  axial  force  coefficient  is  assumed  to  be  contributed  entirely  by  the  body. 

The  total  normal  force  coefficient  is  given  by 

(C^n)tot  ~  sin  0 

+  ICnt2  cos 0|  S.j./S^^j.  +  Ig.j.  (60) 
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C„T2 


Vjj  sin-i 


Figure  12.  Roll  Angle  and  Fin  Load  Definitions 
(Looking  Forward) 


where  S-j-  is  the  fin  planform  area,  S,gf  is  the  body  reference  area,  is  the  body-alone  normal 
force.  Coefficients,  are  the  individual  fin  loads  in  the  presence  of  the  body  and  other  fins, 

and  IgT  is  the  tail-to-body  carryover  normal  force  coefficient. 

The  longitudinal  center  of  pressure  (from  the  nose)  for  the  entire  configuration  is  given  by 

l^cp^TOT  ~  {^cpB  ^NB  ^  y  j  +  ^^cp^T3  ^NT3  ^ 

+  COS0l(x,p)-r2CN^2  +  (x,p)T4  CnT41]  Sr/Sef 
lBrXcpl}ACN)-|-OT 

All  Xj.p  values  are  given  in  calibers;  x^pp  is  the  body-alone  x^p,  (X(-p)Ti  are  the  individual  tail  x^p 
values  in  the  presence  of  the  body,  and  x<,p|  is  the  center  of  pressure  for  the  carryover  load. 

If  one  neglects  the  small  body  frictional  roll  coefficient  and  the  tail  carryover  onto  the  body, 
one  can  estimate  the  roll  coefficient  as 

Cs  =  ^  ‘^nt.  sin  [(i2  -  i  +  1)  4]  n 

i=  1  L  “  J  L  '  /  /J 

Tail  carryover  data  were  unavailable  for  side  force  and  side  moment  coefficient  computa¬ 
tions.  Details  of  the  individual  component  functional  forms  are  available  in  Reference  35  and  will 
not  be  presented  here. 


IMPROVED  DYNAMIC  DERIVATIVE  PREDICTION 


Ml  ■  IW.III  I . . 


A  portion  of  this  work  was  performed  under  contract  by  L.  Ericsson  of  Lockheed  Missile 
and  Space  Co.  (LMSC).^®"'**^  Only  C^^^  +  Cn,a  prediction  is  considered.  In  fact,  the  Cp^a  or  con¬ 
tribution  term  of  constant  vertical  acceleration  is  not  considered.  Restrictions  on  the  LMSC  code 
are 

1 .  Body  or  body-tail  configurations 

2.  Aspect  ratio,  ,  of  the  wing-alone  limited  to  less  than  2.3 

3.  Initial  cylindrical  radius  for  a  spherically  blunted  or  truncated  body  limited  to  less  than  0.25 

calibers  _ 

4.  Contribution  of  afterbody  and  boattail  or  flare  are  neglected  at  Moo  ^  M*  (see  Equation  64) 

5.  Dynamic  derivative  evaluated  in  the  neighborhood  of  a  =  0. 


At  low  Mach  numbers,  one  has  the  option  of  using  the  older  method.  At  higher  Mach  numbers, 
the  contribution  of  the  lifting  surfaces  needs  to  be  supplemented  in  regions  where  the  LMSC 
routine  is  invalid  or  not  used. 


LMSC  Body  Pitch  Damping 

In  subsonic  flow  C^^  +  C^a  is  given  by  a  relation  based  on  slender  body  theory 

(Cmq  +  C^a)B  =  -4(0.77  +  0.23  H2)(eV0.77  +  0.23  M2,  -  (63) 


Here  C  is  the  body  length  in  calibers  and  x^g  is  the  moment  and  rotation  reference  point  from 
the  nose  in  calibers.  The  range  of  validity  is  M.^  <  1 . 

Hypersonic  flow  approximations  are  applied  above  a  certain  Mach  number,  M^  =  M*  associated 
with  an  effective  hypersonic  similarity  parameter  of  0.4  such  that 


M*  =  0.4  CSC  0* 


(64) 


e* 


tan“i 


(65) 


where  Lj^  is  the  length  of  the  nose  in  calibers.  If  M*<  1 .5,  then  M*=  1 .5  is  chosen.  For  Moo>  M*, 
an  embedded  Newtonian  approximation  is  used.  For  1  <  M„  <  M*  a  linear  interpolation  is  used 
where 


M*-M„ 


M*  -  1 


/  M.^  -  1  \ 

_  j)  ^md)B(M„,=M*) 


(66) 


For  M„  >  M*  the  computation  proceeds  as  follows: 


CpoC-y 


(x-x.„  +  RR')2rR' 


(1  +R'2) 


where  R  and  R'  are  the  body  radius  in  calibers  and  x  derivative,  respectively,  x  is  in  calibers,  6  is 
the  local  surface  angle  associated  with  R'  (see  Figure  5),  and  Cp^  is  the  stagnation  pressure  behind  a 
bow  shock 


C  =  — 

po  y 


±2  f,  +_L5_  J_' 
+  1  L  ^  Mi. 


Equation  68  is  approximate  for  large  Mach  numbers.  Also 


1.01  +  1.31  [log  (10  Moo  sin  6)1  M,„  sin  6  >  0.4 


M,^  sin  6  <  0.4 


Note  that  flared  bodies  are  not  considered. 


LMSC  Wing-Body  Pitch  Damping 


For  Moo  ^  1  >  slender  wing-body  theory  is  used  where 

(^mq  ^md  )w  ~  (^mq  ^md)wB  (^mq  Cmd^B  (70) 

where  the  subscript  B  refers  to  the  body,  w  to  the  wing,  and  wB  to  the  wing-body  combination. 
Finally 


«7n,q  +C„,d)wB  =  -4  -pr  (1  -0.23/J^)  \  \x^  v/1  -  0.23^'  -  x 


-  x,g  -  C/D  1 


Here  D  is  the  reference  diameter  of  the  body,  D,,,  is  the  mean  diameter  of  the  body  near  the 
wing,  is  the  root  chord,  and  is  the  x  distance  in  calibers  from  the  nose  to  the  beginning  of 
the  root  chord.  Other  variables  are 

P  =  Vn  -  Ml\  (72) 


C,,,(l  +\) 


=  aspect  ratio  = 


(73) 


X  =  =  taper  ratio 


(74) 


C/D  =  C,^./D 
1  +  2X 


-[ 


1  -  (1  -X)  - 


bw 

2C„ 


tan  n 


V  = 


x,^/D  = 


3(1  +X) 
x^D  +  C 


- !Z  _  f /D 

D  D 


Kn.a  =  2y 


1  +  (A/3/4)  +  ,/l  + 


(¥) 


(75) 

(76) 

(77) 

(78) 


and 


A  =  (A/A,j,)A^ 


The  parameter,  A,  is  the  aspect  ratio  for  the  planform  obtained  by  extending  the  leading 
and  trailing  edges  to  the  body  centerline,  i.e. 


(A/A^.)  = 


1  +  D^IK 


I  + 


1  ~  X 
bj.  1  +  X 


For  supersonic  flow,  1  <  <  V'  +  (A.,j,/4)"^  Fquation  70  again  applies 


+C^a)wB  =  -4(b,/D)2K 


VD-x,,  -C/d(i-v1^ 


1  - 


'XwD  + 
D 


v2-l 


-  X 


eg 


^cg 


<e„,  (C„.)b 


For  hypersonic  flow,  >  Vl  +(Aw/4)^ 


(80) 


(81) 


(C^q+C,d)w  =  %  J 


1  +  X 


+  4^*(f*-a*) 


+  1  a*2  -  (1 -X)(2f*2-8/3{*  a*  +  a*2) 


k*  _  Cfv 


r  = 


1  ^rw 


D 


’‘eg  2  D 


(82) 


(83) 
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f*  is  the  dynamic  pressure  ratio  across  the  bow  shock  and  is  defined  as 


"1.0 

Kn  <  1.25 

Kn-  1.2S 

f*  =  < 

lo" 

12.669724 

1.25  <  Kn  <  11 

(85) 

..17 

A 

z 

Km  =  sin  (86) 

divi  =  sin  '  f'DN/(^7C'po)  (87) 

6f^  is  an  effective  cone  angle  for  the  bow  shock  in  which  the  wing  is  immersed  and  is  the 
nose  wave  drag  where 


Cdn 


C^Cpo  sin^SRR'  dx 


(88) 


Supplemental  Dynamic  Derivative  Prediction 

The  LMSC  work  was  assumed  to  be  able  to  replace  the  empirical  prediction  of  Reference  1 7 
for  the  body-alone  pitch  damping  coefficient.  Since  the  LMSC  prediction  neglects  the  effects  of 
an  afterbody  and  a  boattail  at  supersonic  Mach  numbers,  the  method  of  Reference  17  often  pre¬ 
dicts  better  values  for  the  pitch  damping  (as  will  be  shown  later).  For  this  reason  and  because  of 
the  other  limitations  indicated  in  the  last  two  sections,  the  pitch  damping  LMSC  prediction  for 
body-alone  and  body-tail  is  currently  an  option  in  the  aeroprediction  code.  Even  when  the  LMSC 
package  is  chosen,  the  various  elements  of  the  package  will  be  bypassed  when  any  of  the  code 
restrictions  are  not  met. 

At  higher  Mach  numbers,  where  M^  >  Mg,  the  potential  models  for  determining  the  lifting 
surface  roll  damping  and  pitch  damping  derivatives  are  assumed  to  be  invalid  and  are  supplemented 
by  the  strip  theory  methods. 

For  a  wing,  the  strip  loading  for  Cgp  and  C^q  is  assumed  to  be  proportional  to  the  local 
chord  such  that 


_ N(CNa)w _ 

D=(^)  %=  (l*M. 


C(y)  dy 


N(C|^a),y _ 

6D^b^(l  +X) 


(b^  +D^)^X  -  (D^)^ 


1  (1  -X) 
4  b. 


[(b^  +DJ'*  -  Dtl 


} 


(89) 
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For  the  wing-alone  case,  D  is  replaced  by  for  determining  Cgp .  Equation  89  predicts  high 
values  of  Cjp .  This  expression  was  subsequently  replaced  by  an  expression  that  assumes  the  loading 
goes  to  zero  at  the  tip  as  Vl-Ty/b^^T^  0  such  that 


Cep 


N(CNa)w 


+  u 


(90) 


where  is  the  reference  length  and  N  is  the  number  of  fins. 


,  _  jr  (IX) 

*>  '  4  3 


1  (1-X)7r 

‘2  “  3  “  16 


,  _  JL  ^(1 

It  “  t  ^ 


3  16  15 

Similarly,  for  zero  roll  angle  and  zero  angle-of-attack  position 


4(CNa)v 


^^TTTTTT;  -  y  tan  A,  -  X)"  dxdy 

-)Crw(l  +  i)  •'o 


•  mq 


D2(b^/ 


J _ (CnqiIw _  J  ] 

3  (h,^/2)C,,,(l  +  X)D^|tan  A, 


(Xc„  -  x*)^ 


~(xcp  “  tan  A,  j 


1 


tan  A,  - 


Cfw  ( 1  ~  X) 

K/2 


j(Xcp  Cf,j,) 


(Xj.p  ^Cf\ 


-  b^/2  tan 


(91) 


Here  (C^alw  is  the  wing-alone  Cn®  for  two  fins  and  A,  =  Ale  **  the  leading  edge  sweep  angle. 
For  the  wing  alone  case,  D  is  replaced  by  C„f  =  2/3  +  C,^).  Equa¬ 

tions  90  and  9)  are  applied  to  both  the  canard  and  the  tail  and  added  together.  Interference  is 
neglected  as  was  the  case  with  the  high  Mach  prediction  of  F.  De  Jamette. 

For  the  wing-alone  potential  problem  for  constant  roll  rate,  the  loading  is  zero  at  the  root  and 
the  tip.  However,  for  strip  theory,  the  lift  problem  loading  does  not  vanish  at  the  wing  tip. 
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The  Cjip  and  models  given  in  Equations  90  and  91  will  be  evaluated  in  more  detail  later. 

IMPROVED  TRANSONIC  NOSE  WAVE  DRAG  PREDICTION 
Background 

In  Reference  1  th“  inviscid  axial  force  was  assumed  to  be  contributed  by  a  linear  superposition 
of  component  parts  of  the  body  in  the  transonic  flow  regime.  The  nose  wave  drag  was  based  on  a 
fit  to  a  small  disturbance  potential  model  solution  for  a  family  of  pointed  tangent  ogive  noses. 
The  nose  wave  drag  was  fitted  as  a  function  of  Mach  number  and  nose  length.  The  effect  of  nose 
bluntness  was  ignored. 

A  separate  nose  wave  drag  uncoupled  from  the  boattail  and  wake  effects  is  valid  for  supersonic 
Mach  numbers.  It  is  also  valid  for  subsonic  Mach  numbers  where  a  long  cylindrical  afterbody 
follows  the  nose  section. 

When  the  cylindrical  afterbody  is  not  long,  the  boattail  wake  and  nose  flow  fields  are  strongly 
coupled. 

Currently,  the  boattail  wave  drag  is  based  upon  a  small  disturbance  potential  solution  that 
assumes  a  long  cylindrical  afterbody  ahead  of  the  boattail.  and  a  supersonic  Mach  number  close  to 
1.0.  The  boattail  drag  is  assumed  to  decay  linearly  to  zero  at  M^^  =  0.95  from  the  predicted  value 
at  M«  =  1,05. 

The  base  drag  prediction  is  empirical  and  is  determined  by  the  expression  C^g  =  ~CpBA(Moo) 
(RB/Rj^.f)^  where  Cpg^(Moo)is  the  base  pressure  coefficient  for  a  long  afterbody  with  no  boattail. 
This  empirical  relation  is  approximate  and  is  valid  for  very  small  boattail  angles  and  a  long  afterbody. 

For  a  body  with  strongly  coupled  nose,  boattail,  and  wake  effects,  it  would  probably  be  best 
to  vary  the  sum  of  the  wave  drag  and  base  drag  for  Mach  numbers  below  M^  =  1 .0  in  a  manner 
similar  to  the  variation  of  CpgAfMoo).  However,  this  approach  was  not  implemented  due  to  a  lack 
of  data. 

The  remainder  of  the  discussion  deals  with  the  case  of  a  separable  nose  wave  drag. 

In  Reference  41,  it  is  shown  that  the  variation  of  the  nose  wave  drag  is  shape  dependent.  For 
noses  with  a  finite  angle  at  the  shoulder,  the  pressure  drag  does  not  go  to  zero  as  the  Mach  number 
is  decreased  due  to  the  effect  of  the  flow  separation  at  the  shoulder  on  the  pressure  distribution 
on  the  nose. 

In  Reference  41,  it  is  indicated  that  the  nose  shapes  of  interest  are  bounded  approximately 
by  families  of  cones  and  ogives.  This  is  the  approach  taken  here. 
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Wave  Drag  Prediction  for  a  Family  of  Blunted  Tangent  Ogives 


This  work  was  performed  under  contract  by  Nielsen  Engineering  and  Research  (NEAR). 
Initially  the  wave  drag  was  to  be  based  on  interpolation  between  tabulated  wave  drag  values  for  a 
family  of  blunted  tangent  ogives.  The  matrix  of  independent  variables  for  this  situation  is 

=  0.7,  O.9.S.  l.O.S.  1.: 

(nose  length)  =  0.75,  1.0,  l.:.S,  1.5.  2.0,  3.0.  4.0.  5.0 

(nose  radius)  R^  =  0.  0.05,  0.1. 0.15.  0.2,  0.25,  0.30,  35,  0.4.  0.45,  0.5 

The  wave  drag  is  based  on  pressures  generated  by  a  time  asymptotic  solution  of  the  Euler  equations. 
Details  are  reported  in  Relerence  42.  Actual  computations  were  made  at  R^.  =  0.025.  0.05,  0.15. 
0.25,  0.35.  and  0.5;  and  =  l.f .  3.0.  «.0.  Ihe  computations  at  Rnj  =  0.5  are  of  course  independ¬ 
ent  of  nose  length.  The  lull  geometric  matrix  was  determined  by  cross  plotting  and  graphical 
extrapolation  and  interpolation.  Computations  for  short  noses  with  small  bluntness  did  not  agree 
with  the  full  potential  model  of  Reference  43  or  data.  Therefore,  portions  of  the  wave  drag  matrix 
were  recomputed  using  the  Soutli-Jameson  code.'^^  The  new  matrix  generated  is  previously 
unpublished  data  and  is  given  in  Table  I  for  2R,^,  and  2Lj^ . 

Eirst.  for  vanable  geometric  parameters,  the  wave  drags  at  a  given  Mach  number  are  deter¬ 
mined  by  a  four-point  Lagrange  interpolation.  Next,  for  variable  nose  bluntness,  the  wave  drags  at 
a  given  nose  length  are  determined  I'y  a  five-point  Lagrange  interpolation.  Finally,  fora  given  nose 
bluntness,  the  final  wave  drag  is  computed  by  a  five-point  Lagrange  interpolation. 

Two  problems  arise  m  implementing  this  wave  drag  algorithm.  For  blunter  noses,  the  wave 
drag  is  still  significant  for  Mach  numbers  below  M.„  =  O.K.  In  this  case,  a  quadratic  decay  to  zero 
wave  drag  at  M^  =  0.5  is  assumed.  For  pointed  ogives,  the  apex  angle  may  lead  to  a  detached  shock 
at  M„  =  1.2  for  Lf,;  <  3.  .Small  ilisturbaiice  potential  solutions,  however,  are  possible  between 
M<„  =  1.2  and  the  detachment  Mach  number  provideil  Ln,  is  limited  to  values  greater  than  about 
1.75.  For  Lm  =  1.75.  the  initial  angle  is  approximately  equal  to  the  Mach  wave  angle  at  detachment. 
This  coincides  with  a  Macli  number  of  the  order  1.5.  !  igvire  13  ^hows  the  approximate  variation 
with  Mach  number  tor  a  few  pointed  tangent  ogives.  I  or  M.^  <  1.2.  C|,  is  given  by  the  NEAR 
tabulated  data.  For  M^  >  1.5  computations  are  established  using  a  second-order  shock-expansion 
computation  program  described  in  an  earlier  section.  Note  that  the  second-order  shock-expansion 
method  may  not  be  applied  below  conic  shock  attachment  Mach  numbers  with  any  accuracy.  A 
potential  solution  of  the  type  described  in  Relerence  I  will  fill  the  gap  between  M^  =  1.2  and  the 
attachment  Mach  number.  For  L,,^  <  1  J5.  however,  the  initial  slope  will  be  greater  than  the  Mach 
wave  angle,  somewhere  between  M.„  =  1.2  and  the  sliock  attachment  Mach  number. 


Table  1 

•  1 1)  = 

KRn  ' 

)  for  Blunted  Tangent  Ogives 

2Rn 

2Ln 

0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

= 

^  0.8 

1.5 

0.040 

0.038 

0.037 

0.035 

0.032 

0.028 

0.023 

0.018 

0.015 

0.030 

2.0 

0.025 

0.022 

0,016 

0.010 

0.005 

0.002 

0.000 

0.000 

0.000 

0.012 

2.5 

0.010 

0.006 

0.004 

0.002 

0.001 

0.000 

0.000 

0.000 

0.002 

0.015 

3.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.001 

0.005 

0.019 

4.0 

0,000 

O.UOO 

0.000 

0.000 

0.000 

0.000 

0.002 

0.007 

0.012 

0.028 

6.0 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.002 

0.008 

0.023 

0.045 

8.0 

0.000 

0.000 

0.000 

0.002 

0.005 

0.010 

0.015 

0.022 

0.037 

0.055 

10.0 

0.000 

0.000 

0.000 

0.002 

0.005 

0.010 

0.010 

0.024 

0.042 

0.062 

= 

0.95 

1.5 

0.160 

0.197 

0.203 

0.203 

0.198 

0.189 

0.180 

0.167 

0.154 

0.150 

2.0 

0.080 

0.137 

0.145 

0.138 

0.130 

0.118 

0.103 

0.085 

0.070 

0.122 

2.5 

0.040 

0.080 

0.060 

0.086 

0.075 

0.070 

0.063 

0.047 

0.038 

0.1  10 

3.0 

0.015 

0.052 

0.062 

0.062 

0,050 

0.042 

0.034 

0.028 

0.033 

0.1 16 

4.0 

0.000 

0.15 

0.027 

0.027 

0.020 

0.015 

0.012 

0.008 

0.050 

0.132 

6.0 

0.000 

0.000 

0.005 

0.002 

0.002 

0.003 

0.010 

0.027 

0.075 

0.152 

8.0 

0.000 

0.000 

0.002 

0.002 

0.004 

0.015 

0.018 

0.055 

0.088 

0.160 

10.0 

0.000 

0.000 

0.000 

0.002 

0.009 

0.015 

0.025 

0.055 

0.090 

0.164 

Moc  -• 

1.05 

1.5 

0.280 

0.333 

0.335 

0.333 

0.326 

0.322 

0.313 

0.300 

0.285 

0.270 

2.0 

0.200 

0.268 

0.272 

0.268 

0.262 

0.252 

0.238 

0.221 

0.207 

0.249 

2.5 

0,155 

0.220 

0.225 

0.222 

0.213 

0.201 

0,185 

0.170 

0.162 

0,248 

3.0 

0.135 

0.165 

0.175 

0.175 

0.167 

0.152 

0.135 

0.130 

0.157 

0.249 

4.0 

0.1  10 

0,130 

0.133 

0.130 

0.122 

0.1 1  1 

0.100 

0.102 

0.161 

0.252 

6.0 

0.078 

0.080 

0.076 

0.070 

0.064 

0.062 

0,070 

0,105 

0.170 

0.257 

8.0 

0.055 

0.050 

0.045 

0.043 

0.047 

0.055 

0.070 

0,105 
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Figure  13.  Pointed  Tangent  Ogive  Wave  Drag  Coefficients  (a  =  0) 


Wave  Drag  Prediction  for  a  Family  of  Blunted  Cones 

As  indicated  in  the  introduction  and  at  the  beginning  of  this  section,  the  nose  wave  drag  is 
shape  dependent.  The  approach  taken  here  is  that  the  nose  is  either  a  blunted  tangent  ogive  or  a 
blunted  cone.  The  only  alternative  would  be  to  use  a  code  such  as  that  presented  in  Reference  43 
to  generate  families  of  solutions  for  various  nose  shapes. 

Figure  14  shows  the  pressure  drag  versus  Mach  number  variation  for  various  pointed  cones. 
The  pressure  drags  are  established  by  integration  of  the  pressure  data  from  Reference  45.  Note 
that  the  transonic  flow  regime  or  shock  attachment  does  not  start  until  Mj„  >  1 .2  for  angles  above  15°. 

Data  were  available  up  to  6  =  40°.  However.  6  =  25°  is  about  the  maximum  practical  angle 
for  many  applications.  Note  that  the  pressure  drags  are  higher  than  those  for  tangent  ogives  of  the 
same  length.  Noses  such  as  secant  ogives  may  be  fitted  by  a  cone  with  the  effective  cone  angle 
approximately  et|ual  to  the  mean  value. 

The  total  forebody  drag  (pressure  drag  plus  viscous  drag)  is  relatively  independent  of  blunting 
for  Mach  numbers  below  shock  attachment.^*’  Data  in  Reference  46  are  for  cones  instead  of  cone- 
cylinders.  The  drags  are  different  for  the  two  cases;  however,  it  is  assumed  that  the  effect  of 
blunting  is  the  same. 


1  or  Mach  numbers  below  Mo,,  =  0.‘).  separation  effects  at  the  cone-cylinder  junction  affect  the 
pressure  on  the  nose.  As  the  Mach  number  is  lowered,  the  pressure  drag  does  not  go  to  zero.  The 


Figure  14.  Transonic  Pressure  Drag  for  Pointed  Cone-Cylinders 


residual  drag  was  termed  viscous  drag  in  Reference  1.  However,  the  magnitude  indicated  in 
Reference  I .  is  incorrect,  since  the  data  used  are  based  on  cones  rather  than  cone-cylinders. 

The  drag  behavior  below  Moo  ~  0-7  was  established  by  the  assumed  asymptote 

Ap  v  Apci  ^  '-  Apl'“ 

All  of  the  undetermined  coefficients  may  be  determined  from  C^p  data  at  Moo  =  0.7.  0.8.  and  0.9. 
Table  2  gives  the  resulting  fit  parameters. 

f-'igure  15  is  a  plot  of  C,^p„  against  5  compared  with  the  asymptote  from  Reference  1.  which 
is  assumed  to  apply  below  Moo  =0.8  only.  The  asymptotic  behavior  determined  here  is  indicated  in 
Figure  14  by  dashed  curves.  For  all  practical  purposes,  the  asymptote  for  C.  is  reached  by  Moo  =  0.5. 

A  p 


Table  2,  l  it  I’arameters 


6 

^  Apti 

^  Apl 

A 

5 

0.00041 
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15 
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Figure  15.  Separation  Pressure  Drag  on  a  Cone-Cylinder 


The  pressure  drag  for  cones  may  then  be  given  by  table  lookup  between  the  limits 
0,5  <  M„  <  1.5.  By  limiting  the  angle  to  20°.  the  Mach  number  range  becomes  0.5  <  <  1.2. 

Method  of  Analysis  fora  General  Nose  Shape 


At  the  beginning  of  this  section,  it  was  indicated  that  the  wave  drag  is  dependent  on  shape. 
It  remains  to  estimate  the  wave  drags  for  nose  shapes  that  are  neither  blunted  tangent  ogives  or 
blunted  cones.  Unfortunately,  not  enough  numerical  or  test  data  are  available  to  test  any  empir¬ 
icism.  It  would  be  advantageous  to  consider  a  power  law  base  as  well  if  these  data  were  available. 

One  possibility  is  a  linear  interpolation  as  given  in  the  relation 

^  KCaw)T.O.(R'l  -K)  + 

R;  .  R',  "  ^ 

where  (CJ^^,)^^  o.  the  tangent  ogive  value  of  wave  drag,  6,  is  the  shoulder  angle,  and  6]  is  the 
initial  nose  angle  for  a  pointed  or  truncated  nose.  Here  Sj  is  the  initial  angle  at  the  end  of  the 
spherical  cap  for  a  spherically  blunted  nose  and  (C^,j. )^.(6s)  is  the  cone  wave  drag  evaluated  at  6^. 
R'  and  R'  are  slopes  corresponding  to  6|  and  6^. 

Thus,  when  the  shoulder  angle  is  zero,  a  tangent  ogive  value  is  computed.  When  the  shoulder 
and  initial  angle  are  equal,  the  cone  value  is  computed.  For  a  truncated  body,  the  radius  R^  is 
assumed  to  be  that  of  a  pseudosphere  that  could  be  added. 
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Uiitortimately.  not  onougti  experimental  or  numerical  data  are  available  to  test  the  heuristic 
expression. 

IMPROVED  TRANSONIC  NORMAL  FORCE  PREDICTION 
Background 

In  Reference  I,  and  0,,,^  were  assumed  to  be  contributed  by  a  linear  n!  perposition  of 
component  parts  of  the  body  in  the  transonic  flow  regime.  Nose  properties  were  assumed  to  be 
independent  of  other  parts  of  the  body. 

Boattail  properties  were  assumed  to  be  independent  of  the  nose  and  afterbody.  The  nose- 
alone  C|v.(j  was  based  on  the  cone-alone  normal  force  data  and  x^p  is  given  by  slender  body  theory. 

The  afterbody  normal  force,  C^^,  is  given  as  a  function  of  afterbody  length  and  Mach  number 
independent  of  other  parts  of  the  body.  Both  Cj^^  and  x^,p  are  given  by  the  transonic  small  dis¬ 
turbance  theory  of  Reference  6.  F-or  low  Mach  numbers,  C^^  is  predicted  from  data  correlation. 

For  the  boattail,  C^^  is  based  on  data  correlation;  Xj^.p  is  given  by  slender  body  theory.  A  long 
afterbody  is  assumed. 

This  combination  of  methods  used  in  Reference  1  was  necessitated  by  the  lack  of  a  theoretical 
and  experimental  data  base  at  the  time  Reference  1  was  written. 

The  method  of  Reference  47  was  intended  to  replace  that  of  Reference  1.  The  work  was 
performed  under  contract  by  NEAR, 

Analysis  Method 

The  method  of  Reference  47  solves  the  Euler  equations  by  an  implicit  time  asymptotic 
method  at  a  1°  angle  of  attack.  The  resultant  and  were  to  be  fitted  as  a  function  of  body 
geometric  parameters  and  Mach  number.  Geometric  lengths  were  given  in  caliber  dimensions.  The 
parameter  ranges  to  be  considered  are  taken  from  Reference  47  and  are  shown  in  Figure  16  for 
bodies  with  blunted  tangent  ogive  nose  shapes. 

A  number  of  problems  have  arisen  in  implementing  this  method.  For  a  fairly  coarse  grid  of 
4X  (points  along  the  body)  by  1  2  (points  normal  to  the  body)  by  19  (equally  spaced  circumferential 
planes),  approximately  ,^900  CDC  7600  CPU  seconds  arc  required  per  configuration.  A  typical  grid 
taken  from  Reference  47  is  shown  in  Figure  17.  Needless  to  say,  the  computational  time  is 
prohibitive.  For  each  Mach  number,  19  or  20  configurations  were  considered.  In  some  cases,  a  com¬ 
pletely  converged  solution  was  not  obtained.  In  other  cases,  the  grid  was  too  coarse  (for  longer 
bodies). 

The  functional  fit  fomi  is  a  truncated  Taylor  series  with  free  coefficients  fitted  by  least  squares. 
The  functional  form  is  shown  in  Figure  18.  The  configurations  actually  computed  are  given  in 
Tables  3  through  6.  0  in  F'igure  18  refers  to  reference  conditions,  AR^  =  (Rj,^  -  ,25)/.5;  ALj^ 
=  (L^  -  2.5)/5;  AL^  =  (L^  -  2.0)/5;  and  AOj^  =  (0g-5)/15  are  the  quantities  given  in  Figure  18. 
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Figure  18.  Quadratic  Interpolation  Formula  for  Normal-Force-Curve  Slope 


Table  3.  Cases  Calculated  for  =  0.75  and  =  1 


Geometric  Parameters  Euler  Code  Values  Least-Square  Values 


Case 

Rn 

Ln 

La 

(Degrees) 

Lfi 

da 

^cp 

dC^ 

da 

*cp 

I 

0.25 

2.5 

2.0 

5 

1 

1.542 

-0.4714 

1.5776 

-0.3262 

2 

0.05 

2.5 

2.0 

5 

1 

1.750 

0.04893 

1.7747 

0.2146 

3 

0.50 

2.5 

2.0 

5 

1 

1.195 

-1.925 

1.1821 

-1.9698 

4 

0.25 

1.5 

0.25 

5 

1 

1.403 

0.07262 

1.4097 

-0.1311 

5 

0.25 

5.0 

5.0 

5 

1 

1.478 

-2.280 

1.4771 

-2.2896 

6' 

0.05 

2.5 

0.25 

0 

1 

2.393 

1.328 

2.3513 

1.2106 

7 

0.50 

1.5 

5.0 

0 

1 

1.935 

0.8494 

1.9309 

0.8773 

8 

0.50 

5.0 

0.25 

0 

1 

2.085 

1.154 

2.0847 

1.1991 

9 

0.25 

2.5 

2.0 

0 

1 

2.202 

1.147 

2.2267 

1.0143 

10 

0.05 

5.0 

5.0 

0 

1 

2.295 

2.491 

2.2955 

2.5022 

11 

0.50 

1.5 

5.0 

10 

1 

0.6884 

-9.249 

0.705 1 

-8.9992 

12 

0.50 

5.0 

0.25 

10 

1 

0.2472 

-29.78 

0.2479 

-30.1348 

13 

0.25 

2.5 

2.0 

10 

1 

0.9802 

-3.503 

0.9470 

-3.5864 

15 

0.05 

1.5 

0.25 

10 

1 

0.9875 

-0.5368 

0.9872 

-0.6565 

A3 

0.25 

1.5 

2.0 

5 

1 

1.548 

-0.4303 

1.5073 

-0.5041 

AS 

0.25 

2.5 

5.0 

5 

1 

1.622 

-1.326 

1.6094 

-1.3903 

A8 

0.05 

2.5 

2.0 

10 

1 

1.246 

-1.464 

1.2478 

-1.5871 

AlO 

0.25 

2.5 

0.25 

lO 

1 

0.8119 

-2.269 

0.8249 

-1.8410 

All 

0.10 

1.5 

0.25 

0 

1 

2.124 

0.0821 

2.1458 

0.2236 

Table  4.  Cases  Calculated  for  =  0.90  and  =  1 

Geometric  Parameters 

Euler  Code  Values 

Least-Square  Values 

Case 

dC. 

Rn 

Ln 

La 

(Degrees)  ®  da 

^cp 

da 

*cp 

1 

0.25 

2.5 

2.0 

5 

1  1 .498 

-0.8088 

1.4542 

■1.0939 

2 

0.05 

2.5 

2.0 

5 

1  1.689 

0.0848 

1.6834 

-0.0617 

3 

0.50 

2.5 

2.0 

5 

1  1.305 

-1.973 

1.3115 

-1 .8470 

4 

0.25 

1.5 

0.25 

5 

1  1.479 

-0.03936 

1.4685 

-0.3115 

5 

0.25 

5.0 

5.0 

5 

1  1 ,437 

-2.764 

1.4424 

-2.7339 

6 

0.05 

2.5 

0.25 

5 

1  1 .494 

0.3263 

1.6030 

0.6461 

7 

0.50 

1.5 

5.0 

0 

1  2.124 

1.178 

2.1156 

1.1613 

8 

0.50 

5.0 

0.25 

0 

I  2.008 

0.7239 

2.0057 

0.7358 

9 

0.25 

2.5 

2.0 

0 

1  2.240 

1.113 

2.2653 

1.0701 

10 

0.05 

5.0 

5.0 

0 

1  2.302 

2.433 

2.2981 

2.4400 

11 

0.50 

1.5 

5.0 

10 

1  0.6003 

-12.23 

0.6023 

-12.3097 

12 

0.50 

5.0 

0.25 

10 

1  0.4092 

-18.29 

0.4096 

-18.4090 

13 

0.25 

2.5 

2.0 

10 

1  0.9875 

-3.926 

0.9687 

-4.0208 

15 

0.05 

1.5 

0.25 

10 

1  1 .032 

-0.8020 

1 .0293 

-0.8941 

A3 

0.25 

1.5 

2.0 

5 

1  1.388 

-1.751 

1.4182 

-1.3427 

A5 

0.25 

2.5 

5.0 

5 

1  1.562 

-1.808 

1.5659 

-1.7509 

A8 

0.05 

2.5 

2.0 

10 

1  1.369 

-1.125 

1.3755 

-1.0947 

AlO 

0.25 

2.5 

0.25 

10 

1  0.8812 

-2.247 

0.8938 

’  9900 

Al  1 

0.10 

1.5 

0.25 

0 

1  2.366 

0.8041 

2.3546 

.  3441 

Table  5.  Cases  Calculated  for  =  0.95  and  Ly  =  1 


Geometric  Parameters  Euler  Code  Values  Least-Square  Values 


Case 

Rn 

Ln 

La 

(Degrees) 

Lb 

da 

*cp 

££n 

da 

*cp 

1 

0.25 

2.5 

2.0 

5 

1 

1.641 

-0.5432 

1.5336 

-0.9262 

2 

0.05 

2.5 

2.0 

5 

1 

1.833 

0.2568 

1.8396 

0.1436 

3 

0.50 

2.5 

2.0 

5 

1 

1.361 

-2.192 

1.3652 

-2.1112 

4 

0.25 

1.5 

0.25 

5 

1 

1.352 

-0.7351 

1.3685 

-0.8539 

5 

0.25 

5.0 

5.0 

5 

1 

1.543 

-1.964 

1.5733 

-1.8363 

6 

0.05 

2.5 

0.25 

5 

1 

1.652 

0.3024 

1.6443 

0.3821 

7 

0.50 

1.5 

5.0 

0 

1 

2.225 

1.048 

2.1860 

0.9678 

8 

0.50 

5.0 

0.25 

0 

1 

2.118 

0.6828 

2.0962 

0.6421 

9 

0.25 

2.5 

2.0 

0 

1 

2.193 

1.021 

2.3506 

1.2406 

10 

0.05 

5.0 

5.0 

0 

1 

2.330 

2.506 

2.3057 

2.4914 

11 

0.50 

1.5 

5.0 

10 

1 

0.5024 

-20.17 

0.5382 

-18.5729 

12 

0.50 

5.0 

0.25 

10 

1 

0.7661 

-8.846 

0.7820 

-8.5952 

13 

0.25 

2.5 

2.0 

10 

1 

1.252 

-2.260 

1.0861 

-3.4230 

15 

0.05 

1.5 

0.25 

10 

1 

1.323 

-0.4737 

1.3319 

-0.4645 

A3 

0.25 

1.5 

2.0 

5 

1 

1.442 

-1.717 

1.4924 

-1.3205 

A5 

0.25 

2.5 

5.0 

5 

1 

1.632 

-1.639 

1.6299 

-1.6229 

A8 

0.05 

2.5 

2.0 

10 

1 

1.456 

-0.8952 

1.5107 

-0.6385 

AlO 

0.25 

2.5 

0.25 

10 

1 

0.8763 

-2.764 

0.9269 

-2.2428 

All 

0.10 

1.5 

0.25 

0 

1 

2.471 

0.8912 

2.3989 

0.8074 

Table  6.  Cases  Calculated  for  M„  =  1-2  and  Lg  =  1 

Geometric  Parameters  Euler  Code  Values  Least-Square  Values 

Og  dCfj  dC^i 


Rn 

Ln 

La 

L  — ~ 

(Degrees)  da 

^cp 

da 

*cp 

1 

0,25 

2.5 

2.0 

5 

2.136 

0.1233 

2.1773 

0.6058 

2 

0.05 

2.5 

2.0 

5 

1.989 

0.7118 

1.9209 

0.3301 

3 

0.50 

2.5 

2.0 

5 

2.167 

0.0200 

2.1823 

0.0313 

4 

0.25 

1.5 

0.25 

5 

2.740 

0.6582 

2.7247 

0.7038 

5 

0.25 

5.0 

0.50 

5 

2.334 

1.403 

2.3160 

1.2874 

6' 

0.05 

2.5 

0.25 

0 

2.500 

1.494 

2.5538 

1.5969 

7 

0.50 

1.5 

5.0 

0 

2.584 

0.8771 

2.6007 

0.9326 

8 

0.50 

5.0 

0.25 

0 

2.592 

0.7101 

2.6184 

0.7563 

9 

0.25 

2.5 

2.0 

0 

2.562 

0.8747 

2.4646 

0.5613 

10 

0.10 

5.0 

5.0 

0 

2.180 

-1.888 

2.1965 

-1  7688 

1 1 

0.50 

1.5 

5.0 

10 

1.983 

-1.035 

1.9505 

-1.1661 

12 

0.50 

5.0 

0.25 

10 

1.894 

-0.9524 

1.8917 

-0.9216 

13 

0.25 

2.5 

2.0 

10 

1.833 

-0.6771 

1.9130 

-0.2967 

15 

0.05 

1.5 

0.25 

10 

2.176 

0.7843 

2.1743 

0.8177 

A3 

0.25 

1.5 

2.0 

5 

2.247 

-0.0886 

2.2748 

-0.1955 

A5 

0.25 

2.5 

5.0 

5 

2.265 

-0.0367 

2.2831 

0.0177 

A8 

0.05 

2.5 

2.0 

10 

1.751 

0.0354 

1.7516 

0.1103 

AlO 

0.25 

2.5 

0.25 

10 

2.323 

0.7546 

2.2464 

0.4127 

All 

0.10 

1.5 

0.25 

0 

2.846 

1.153 

2.8508 

1.1936 
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The  Taylor  series  is  used  for  the  nose-alone;  afterbody-alone;  and  0.5-,  1.0-,  and  1.5-caliber 
boattails.  In  the  actual  computations,  only  a  2-caliber  boattail  is  considered.  The  loading  on 
0.5-,  1.0-,  and  1.5-caliber  boattails  is  obtained  by  integrating  over  the  truncated  portion  of  the  full 
2-caliber  boattail. 

The  interpolation  procedure  proceeds  as  follows;  At  =  0.6,  the  nose-alone,  afterbody- 
alone,  and  boattail-alone  and  values  are  obtained  by  the  method  of  Reference  1.  Next,  at 

=  0.75,  0.9,  0.95,  and  1.2,  the  nose-alone;  afterbody-alone;  and  0.5-,  1.0-,  and  1.5-boattaiT 
length  contributions  to  Cfja  and  are  determined  from  the  Taylor  series.  At  each  Mach  number, 
the  boattail-alone  contribution  is  determined  by  a  four-point  Lagrange  interpolation  for  the  boattail 
length  of  interest.  For  a  zero  length  boattail,  and  contributed  by  the  boattail  are  zero. 
Finally,  the  nose-alone,  afterbody-alone,  and  boattail-alone  contributions  are  determined  at  the 
Mach  number  of  interest  by  a  five-point  Lagrange  interpolation  in  the  values  at  M^  =  0.6,  0.75, 
0.9,  0.95,  and  1.2.  Since  two  different  methods  are  used,  a  smooth  merge  would  not  be  expected  at 
Moo  ~  0-6-  At  Moo  ^  1-0  Oie  code  as  currently  written  would  transfer  to  a  potential  code  for  the 
low  supersonic  range.  Again,  a  smooth  merge  would  not  be  expected  for  the  current  new  coding  or 
the  older  coding  at  Moo  ~  1-20. 

At  subsonic  Mach  numbers,  the  hemispherical  end  configuration  leads  to  poor  pressure  dis¬ 
tribution  prediction  on  the  aft  end  of  the  body,  especially  for  -  0°.  Some  guidance  from  flow 
visualization  is  needed  in  choosing  a  free  streamline  wake  boundary  which  may  be  approximated  by 
a  solid  surface.  The  next  best  thing  is  an  extended  smoothly  varying  stinglike  wake  free  streamline. 
The  alternative  is  to  couple  a  boundary  layer  and  transonic  computation  such  as  used  in  Reference 
48.  The  model  of  Reference  48  would  have  to  be  extended  to  include  a  recirculation  model  near  the 
base  or  a  knowledge  of  the  free  streamline  shape. 

Examination  of  dC^/dx  in  Reference  47  (see  Figure  19)  indicates  that  this  model  predicts  a 
boattail  loading  at  0^  =  0.  This  is  questionable,  particularly  for  M^^  =  1.2,  for  the  example  shown  in 
Figure  19.  Unfortunately,  the  0^  =  0  numerical  data  are  included  in  the  Taylor  series  fit.  The 
truncated  portions  of  a  2-caliber  boattail  are  not  the  same  as  those  of  0.5  .  1.0-,  and  1.5-caliber 
boattails.  However,  the  prohibitive  computational  cost  led  to  the  current  approach. 

A  proper  evaluation  of  the  current  method  would  ideally  require  extensive  dC^/dx  distribu¬ 
tion  data  or  Cp  data  that  could  be  used  to  determine  dCnj/dx.  The  majority  of  the  data  available 
are  total  load  data  for  complete  configurations.  However,  the  basic  input  data  to  the  fits  need 
to  be  examined  and  evaluated  as  well. 


Evaluation  of  Current  and  Old  and  Cp,Q  Predictions 

One  limit  that  one  could  examine  is  the  case  of  a  long  constant  radius  afterbody.  In  principle, 
the  nose  and  boattail  should  be  isolated  in  this  case.  Indeed,  this  was  the  assumption  of  Reference  1. 

Nose-Alone  Contribution.  Nose  geometry  considered  could  be  listed  as  (see  Tables  3-6)  in 
Table  7. 


Table  7.  Nose  Geometries 


Rn  *-n 


0.05 

1.5 

0.05 

2.5 

0.05 

5.0 

0.1 

1.5 

0.1 

5.0 

0.25 

1.5 

0.25 

2.5 

0.25 

5.0 

0.5 

1.5,  2.5,  5 

Case(s) 

15 

2,  6',  A8 

10  (not  for  =  1.2) 
All 

10  (only  for  M„  =  1.2) 
4,  A3 

1,9,  13,  A5,  AlO 
5 

(7,11),  (3),  (8,  12) 


The  hemispherical  cases  shown  in  Table  7  are  used  in  the  fits.  In  actual  usage,  a  hemispherical 
cap  is  given  by  =  0.5,  =  0.5. 

Values  of  C,,  and  C  ^  for  the  nose  alone  were  taken  from  unpublished  data  obtained  from 
G.  Klopfer.  Computational  data  for  cases  2,  6',  and  A8  are  tabulated  in  Table  8.  Cases  2,  6  ,  and  A8 
have  different  values  of  and  6^.  However,  the  data  seem  to  indicate  that  the  assumption  of 
nose  contribution  independent  of  afterbody  and  boattail  is  not  too  bad. 

Cases  1,  9,  13,  A5,  and  AlO  also  have  a  common  nose.  The  nose-alone  contribution  is  tabu¬ 
lated  in  Table  9. 

Again,  the  trend  toward  nose-alone  values  independent  of  the  aft  portion  of  the  body  is 
apparent.  An  exception  is  for  cases  4  and  A3.  For  a  nose  as  blunt  as  in  cases  4  and  A3,  the  after¬ 
body  and  boattail  apparently  influence  the  nose.  Case  4  has  a  very  short  =0.25,  as  opposed  to 
=  2  for  case  A3. 

However,  the  numerical  data  base  available  to  take  advantage  of  this  result  is  far  too  small. 
The  values  at  *  are  suspect  since  for  case  10  a  value  of  Cfga  =  3.03,  which  is  much  too  high, 

is  obtained  for  the  nose.  Including  dependence  of  Cfga.  and  C^^  on  and  in  the  curve  fit  for 
the  nose,  distorts  the  contribution  of  these  elements. 


Table  8.  Cn„  and  C™^  Values  for  Nose-Alone  Cases  2,  6  ,  and  A8 


CNa 

M„/Nose- Alone 

Case  2 

Case  6' 

0.75 

2.31 

2.28 

0.9 

2.24 

2.17 

0.95 

2.45 

2.46 

1.2 

2.28 

2.28 

ma 


Case  A8 

Case  2 

Case  6' 

Case  A8 

2.14 

-2.86 

-2.79 

-2.76 

2.33 

-2.87 

-2.75 

-3.15 

2.44 

-3.37 

-3.39 

-3.35 

2.25 

-2.97 

-3.14 

-2.95 
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Table  9.  and  Values  for  Nose-Alone  Cases  1,9,  13,  A5,  and  A 10 


M  /Nose- Alone 


Cases  1, 9,  13,  A5,  AlO 


Cases  1. 9,  13,  A5,  AlO 


2.07,  2.07,  2.08,  2.06,  2.07 
2.12,  2.15,  2.16,  2.13,  2.23 
2.38,  2.30,  2.38,  2.41,  2.40 
2.73,  2.78,  2.72,  2.64,  2.80 


-2.04,  -2.04,  -2.05,  -2.01,  -2.01 
-2.01.-2.03,  -2. 11,  -2.00,  -2.25 
-2.58,  -2.29,  -2.58,  -2.57,  -2.52 
-3.21,-3.46,  -3.24,  -3.30.  -3.54 


Table  10.  Cn^  and  C„,a  Values  for  Nose-Alone 
Cases  4  and  A3 


M 


Case  4 

Case  A3 

Case  4 

Case  A3 

0.75 

2.01 

2.05 

-1.43 

-1.48 

0.9 

2.45 

2.72 

-1.95 

-2.1 1 

0.95 

2.47 

2.74 

-2.01 

-2.15 

1.2 

2.94 

2.74 

-2.52 

-2.28 

A  better  approach  would  have  been  to  use  a  potential  model  such  as  developed  in  Reference  49 
for  noses  plus  very  long  afterbodies.  A  possible  fit  would  be 

G(L*,  Or.) 

(CNa)N  =  F(Rn,  Ln)  +  -  ^  (93) 

The  last  term  would  provide  correction  terms  for  short  bodies.  For  supersonic  Mach  numbers 
the  second  term  in  Equation  93  would  not  be  needed. 

Another  question  that  arises  is  to  what  extent  the  nose  loading  is  dependent  on  the  shape  of 
the  nose.  Reference  1  bases  the  nose  loading  on  the  cone-alone  data  from  Reference  46  where 

(FNa)N  =  C,(M^)tan6*  +  (94) 

where  6*  is  the  angle  at  the  shoulder.  Thus,  for  the  family  of  blunted  tangent  ogives,  the  prediction 
's  Na^N  ~  F2(Moo)  or  independent  of  the  nose  parameter  For  tan  6*  -*■  0,  {CNa)^  -*  2.0  is  the 
trend  of  the  data  for  cones.  However,  C2(Moo)  is  above  2.0  for  higher  Mach  numbers;  this  is  just  a 
property  of  the  force  fit. 

Computations  from  Reference  47  and  data  from  Reference  49  indicate  that  (C^alN  for 
blunted  tangent  ogives  is  above  2.0.  Data  from  References  45  and  46  indicate  that  (C^a^N 
than  2.  Reference  46  indicates  that  bluntness  does  not  greatly  affect  (Cfg^^jN  in  the  transonic  flow 
regime  for  cones  for  Mach  numbers  below  the  shock  attachment  Mach  number  for  a  pointed  cone. 
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l  ifniR'  20  Is  a  i.-ompaiisoii  Ih  Iui  cii  iIk-  iniciJialcil  ilala  ol'  Reference  45  for  cone-cy  linders  and 
tlic  values  obtained  in  Releienee  I  loi  cones  alone.  1  he  lit  from  Reference  I  is  low  for  for 
the  three  cone  ancles  plotted  Ihnsever.  I  ..,  preilietetl  by  slender-body  theory  is  too  hijih  com- 
paretl  with  ilata,  except  below  10  cone  ancle.  IVrcentacewise  the  deviation  is  not  too  jireat 
for  x^p/Ls;  .  The  intecrated  data  liom  Releienee  45  has  some  scatter  but  follows  a  general  trend  with 
the  slender-botly  limit  clearly  bcmc  appnKKlied  as  .5  »  0. 

Afterbody-Alone  Contribution  lleie  the  data  trends  and  computational  trends  are  even  less 
clear.  Reference  1  assumes  that  the  carry  over  from  the  nose  rrr  boattail  is  negligible  and  (( 

(x^,p);^/L^  are  a  function  of  afterbody  length  anti  Mach  number  alone.  The  data  of  Reference  1 
also  indicate  that  (f  is  monotomc.  increasing  to  an  asymptote  as  grows.  The  (C’NaU  from 

Reference  I  appears  to  be  representative  of  tile  buddup  on  the  afterbody  for  a  long  afterbody 
with  a  slender  nose  and  a  very  small  ancle  at  the  shoulder. 

For  cone-cylinders,  buildup  is  available  from  Reference  45.  For  ogive-cylinders,  only 

as  a  function  ot  atterbotly  length  (Reference  50)  is  available.  Unfortunately,  nose-alone  data 
are  not  available  in  References  50  or  5  1 ,  which  also  considers  power-law  nose  shapes.  In  Refer¬ 
ences  50  and  51,  the  miiiimutn  afterbody  length  is  4.  Data  from  References  45  and  50  indicate 
that  Cs,'a  continues  to  increase  with  altcrbody  length  in  an  almost  linear  manner  rather  than  going 
to  an  asymptote.  The  magnituile  of  the  linear  slope  cannot  be  accounted  for  by  a  viscous  crosstlow 
correction,  which  is  also  linear  with  l.^.  since  it  would  require  angles  of  attack  of  the  order  of 
10°  or  more,  'fhe  tnagnitude  (('vvtx,  for  a  l(T  cone-cylinder  is  also  larger  in  Reference  45  than  in 
Reference  1.  Fhe  (C  sialv  contribution  lor  a  10"  cone  at  ,M  =  1.1  increases  with  to  a  maxi- 
nuini  at  about  1.5  caliber,  then  down  to  a  minimum  at  about  3  calibers,  and  then  increases  in  a 
steady  linear  manner  with  increasing  1,,^ .  fhe  data  for  Reference  50  are  only  for  of  4  calibers 
or  greater. 

Whereas,  the  available  tiala  aiul  Reference  I  predict  an  incremental  which  is  positive, 

the  computations  of  Reference  47  preilict  both  positive  aiul  negative  values.  The  computations  of 
Reference  47  predict  positive  increments  on  the  afterbody  only  fora  0°  boattail  angle. 

No  comment  w  ill  be  made  on  the  x,  p  prediction  since  the  prediction  accuracy  is  uncertain. 

Boat  tail- Alone  Contribution  Here  the  ilata  available  are  almost  nonexistent.  Reference  1  uses 
a  extremely  limited  data  base  from  essentially  one  source^'  to  establish  an  empirical  prediction. 
Reference  I  indicates  that  the  correlation  is  supplemented  with  data  for  the  5754  improved  and 
175  mm  projectiles.  I  his  wouUl  liias  the  correlation  for  these  projectiles  since  it  is  assumed  that  the 
nose  and  atterbody  predictions  are  correct  .Again,  x^p  is  given  by  a  slender-body  estimate.  The 
expression  given  in  Relerence  1  is  not  correct  ami  should  be 
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Figure  20.  Comparison  of  Cone-Cylinder  and  Cone-Alone  Data  for  a  Pointed  Conical  Nose 


v^tuTc  t  IS  tile  center  ol  pressure  in  calibers,  t  is  the  body  length  in  calibers,  Lg  is  the  boattail 

length  in  calibers,  (V'olig  is  the  boattail  volume.  S^^.,  “  rr  4,  Sg  =  trRg.  and  Rg  is  the  cylindrical 
radius  at  the  eiul  ol  the  body  m  calibers. 

lor  a  linite  angle  boattail,  the  computations  of  Reference  1  predict  a  downloading.  The 
computations  of  Reference  47  predict  positive  and  negative  loading  for  a  0°  boattail,  which  is 
really  a  continued  alterbody. 

Modified  Functional  Form  Fit 

•A  preliminary  analysis  ot  the  presious  paragraphs  indicates  that  both  the  original  and  the 
more  recent  prediction  nietliods  base  some  shortcomings.  The  NFAR  functional  form  is  invalid 
outside  the  tit  range,  namely  5-caliber  afterbody  lengths.  The  older  computation  is  at  least  bounded 
for  longer  afterbodies  Currently,  the  user  has  a  choice  as  to  which  method  to  use  for  prediction 
purposes.  For  afterbodies  longer  than  5  calibers,  the  older  method  is  used. 

Hie  M  AR  code  was  deemed  to  bt  inaccurate  for  supersonic  Mach  numbers  because  of  the 
=  12  estimate  for  the  individual  beidy  components.  The  =  1.2  values  that  are  used  in  the 
final  interpolation  are  currently  determined  from  functional  form  fits  to  computations  generated 
by  an  mii'roved  Van  Dyke  hybrid  potential  model  (Reference  53).  The  new  functional  form  fits 
were  generated  from  a  much  larger  data  base  than  the  original  form  fits.  The  total  computational 
cost  was  about  one  hour  of  (  IK  (i700  execution  time.  The  C'mq  and  x^.p  predictions  at  ~  1- 
compare  well  with  those  ot  the  original  (lotential  computation. 

SUMMARY  OF  NFW  AN ALYShS  MLTHODS 

Ihe  body-alone  analysis  methods  are  siimmari/.ed  in  Figure  21.  The  Mach  region  division 
points  are  norniallv  at  M  ^  =  O  X.  12,  and  somewhere  in  the  region  2  <  Mg  <  3  as  determined  by 
the  usc-r. 

New  methods  for  computing  the  wmg-alone  and  interference  aerodynamics  are  listed  in 
figure  2  2 

finally,  the  new  methods  used  to  compute  the  dynamic  derivatives  are  summarized  in 
figure  23  file  input  option  refers  to  use-  ol  I  M.S(  or  other  prediction  methods  of  0^,^  +  Cmd 
body-alone  or  body-tail  contigurations 

At  high  angles  ol  attack,  the  body-alone  or  beidy-tail  computations  are  based  on  the  empiri¬ 
cism  of  Reference  35  Only  static  derivatives  are  determined. 

INDIVIDUAL  MFTHODS  FVALUATION 

In  this  section  the  individual,  new  analysis  methods  will  be  evaluated  by  comparison  with 
other  analytical  methods,  numerical  methods,  and  experimented  data  where  available.  No  com¬ 
parison  lor  the  older  methods  will  be  considered  except  where  applicable. 
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Figure  21 .  New  Methods  for  Computing 
Static  Body-Alone  Aerodynamics 
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Figure  22.  New  Methods  for  Computing  Wing-Alone 
and  Interference  Aerodynamics 
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Figure  23.  New  Methods  for  Computing  Dynamic  Derivatives 


HI(,H  M\(  H  NUMBLR  STATIC  ALROUYNAMK  S 


One  consuii.'i;iti()ii  is  the  delermiiiaiioii  of  tlie  Macli  luiiiiher  M  =  M^,.  helow  which  potential 
theory  will  be  used  and  above  which  the  seeoiul-onler  shock-expansion  and  strip  theories  will 
be  used. 

Budy-.Alone  Inviscid  Aerocly  namics 

I  he  High  Mach  body-alone  routine  was  evalualeil  and  compared  with  the  potential  codes  from 
References  1  and  st  and  with  experimental  data.  Space  litnitations  permit  only  a  small  portion  of 
tliis  evaluation  to  be  presented  here.  Originally  it  was  hoped  that  the  High  Mach  body-alone  code 
could  be  extended  down  tt)  M^  =  1.2.  i.e.,  to  be  used  throughout  the  supersonic  How  regime.  This 
did  not  prove  to  be  feasible. 

Although  purely  inviscid  data  lintcgrated  pressure  ilatatarenot  readily  available,  data  for  two 
contigurations.  a  blunt-nosed,  flared  body  and  a  ll.-X  blunted  conical  nose,  for  which  integrated 
pressure  data  arc  available  were  taken  from  Reference  28.  (  omparisons  made  between  these  data, 
the  results  from  the  High  .Mach  code,  and  the  two  potential  codes  of  References  1  and  53  are  given 
in  f  igures  24  and  2.5.  As  can  be  seen,  the  High  .Mach  code  does  a  creditable  job  for  all  functions 
for  these  particular  configurations  and  Mach  number  range  shown.  In  general,  the  drag,  .  pre¬ 
diction  was  found  to  show  belter  agreement  than  that  achieved  for  the  and  x^.p  predictions. 

f  rom  Reference  44.  the  forehody  drag,  .  was  available  for  blunted  tangent  ogive-cylinder 
bodies  with  afterbody  lengths  of  6  cylinders.  Comparisons  of  (\|  for  four  noses  are  given  in 
f  igure  2().  Note  that  the  High  Mach  prediction  accuracy  for  (\|.  is  quite  poor,  below  M^  -  1.5 
.md  improves  with  Mach  number,  for  the  low  supersonic  Mach  range,  the  methoi^  of  Refer^ce  1 
is  best  for  blunted  tangent-ogive  noses.  The  normal  force  derivative.  comparisons  tor  the 

same  four  bodies  are  presented  in  f  igure  27,  while  x^.p  downstream  from  the  nose-afterbody  junc¬ 
tion  is  given  in  figure  28.  fhese  comparisons  indicate  that  the  prediction  using  the  High  Mach 
methodology  degrades  sharply  with  bluntness:  none  of  the  methods  appear  to  compare  well  with  the 
data  in  predicting  x^.p.  Also,  the  trend  for  x^.p  as  predictcl  by  the  High  Mach  code  is  incorrect  for 
the  blunter  configurations  fhe  poorest  prediction  is  for  highly  curved  bodies  and  bodies  with 
negative  slo|ies. 

One  clear  conclusion  ic  that  the  High  Mach  body-alone  prediction  should  not  be  used  below 
M,^  2.  However,  no  othei  rapid  computational  choice  currently  isavailabic  for  predictingC^^  and 

x^.p  at  high  Mach  numbers.  In  the  code  the  dividing  point  between  high  Mach  number  and  low 
supersonic  Mach  number  is  an  input  choice  where  a  range  of  2  .M^'  <  3  is  recommended. 

VVing-.Alone  Inviscid  Aerndynanmw 

for  the  range  oi  .ippheabihty .  Ihe  strip  method  predicts  almost  constant  |3(\  and  jSCfga 
given  or  thickness-lo-chord  r.ilio.  lo  properly  evaluate  the  method,  more  exact  solutions 

or  experimental  data  would  be  leiiuired,  t  Inlorlunately.  neither  was  available.  Instead,  comparisons 
with  the  potential  solution  of  Reference  2  were  made  for  M^  5*  2  as  showni  in  Figure  24  for  j3C;vjQ. 
Ihe  cominirisons  are  shown  lo  improve  with  increasing  aspect  ratio  and  Mach  number,  as  would  be 
ex|iected.  figure  3(1  compare  ^  ‘  IK.  for  ,1  symmeirical  diamond  (double  wedge)  airfoil  where 
fu  ^  r«  for  all  .ases  shown.  Again,  as  one  wouUl  expect,  better  agreement  occurs  for  the 

larger  aspect  ratios.  Although  the  ( \,,  prediction  is  high,  interference  effects  are  not  considered 
anil  hence  in  m.my  cases  the  pre  liclion  o|  at  high  Mach  numbers  is  quite  adequate. 
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11.5°  BLUNTED  CONE  NOSE 
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(c»  Ln  =  1.859.  Rn  =  0.375.  a  =  0 
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(d)  Ln  =  2.357.  Rn  =  0.250.  a  =  0 
Figure  27.  Blunted  Tangent  Ogive-Cylinder.  (Continued) 
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(c)  Ln  =  2.357.  Rn  =  0.250,  a  =  0 
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(d)  Lfi  =  1.859.  Rn  =  0.375.  a  =  0 

Figure  28.  Blunted  Tangent  Ogive-Cylinder,  x^p  (From  Body-Nose  Junction)  (Continued) 


Complete  Configuration  Aerodynamics 


l  igure  31  ctimparcs  dala-^'*  ami  the  High  Mach  computation  for  a  typical  body-tail-canard 
configuration.  1  he  configuration  shown  in  Reference  54  has  body  strakes  not  accounted  for  by 
the  present  combined  theory,  which  partially  accounts  for  the  higher  experimental  values  of  C^^^. 

Data  were  also  available  at  =  4.65  for  the  configuration  sliown  in  Figure  32.  The  body  has 
a  spherical  cap  with  a  0. 1-calibei  radius,  followed  by  a  conical  frustum  3.1]  calibers  long  and 
5  74  '  inclination,  followed  by  a  conical  frustum  2.13  calibers  long  of  inclination  2.36°.  followed  by 
a  constant  radius  afterbody  of  2.67  calibers,  and  ending  with  a  2-caliber.  10°  inclined  Hare.  The 
wing  cross  section  is  a  rectangular  slab  0  t)()2  calibers  thick  with  a  leading  edge  sweep  of  70°  and  a 
trailing  edge  sweep  ot  0°.  I  he  root  cliord  and  the  tip  chord  are  2.46  calibers  and  0.44  calibers, 
respectively.  I  he  nominal  reference  diameter  is  0.25  ft.  Reference  54  indicates  that  the  end  of  the 
body  IS  a  skirt  rather  than  a  tlare.  I  his  could  acccnint  lor  the  difference  in  drag  values. 


HIGH  ANGLh-OF-ATTACK  AERODYNAMICS 

rile  approach  is  purelv  empirical  as  earlier  indicated.  One  configuration  that  lies  in  the  geo¬ 
metric  and  Mach  number  angle-of-attack  envelope  is  the  Air  Slew  demonstrator  configuration. 
(  omparisons  with  experimental  data  are  given  in  Figures  33  through  35.  The  normal  force  coeffi¬ 
cient.  (  \  .  predictedi  adev|uately  and  x^.p  is  predicteil  best  at  higher  angles  of  attack.  However,  the 
(  \  prevlk  turn  is  o|  the  right  order  of  magnituilc  onis . 

\iu)i1icr  comparison  is  made  tor  a  modified  basic  finner  configuration  (Figure  36).  The 
comparison  is  slightly  better  for  x^p  in  this  case  even  though  this  vehicle  differs  from  the  first  only 
sliglith  in  nose  length.  Note  that  the  aerodynamics  for  complete  configurations  are  predicted  to 
a  -  45  oniv.  1  he  method  ot  Reference  56  gives  predictions  to  a  =  180°  but  no  roll  orientation 
dependent  aerodv  mimics  are  computed,  file  current  method  gives  only  the  individual  component 
aerodv  namics  lor  angles  of  .itt.ick  greater  than  45°.  The  one  exception  is  for  the  axial  force  coeffi¬ 
cient  prediction.  I  he  current  method  assumes  that  the  body  contributes  the  majority  of  the  axial 
torce  or  drag. 

DYNAMIC  DERIVATIVE  PREDICTION 

I  he  IM.SC  pitch  damping  prediction  algorithm  was  compared  very  extensively  with  data  and 
a  purely  empirical  (icneral  I  lectric  prediction  program  (under  contract  to  General  Electric  for  the 
evaluation  I,  (  omparisvms  given  in  Reference  57  are  too  numerous  to  be  included  here.  The  con¬ 
clusion  of  R,  Whyte  of  G,F7-‘'^was  that  the  L,SM(  algorithm  was  adequate  for  most  applications. 

I  or  body-alone  cases.  C,„^|  prevliction  was  nullified  as  earlier  indicated.  When  the  value  of 
(  „„i  prevlicted  is  less  negative  than  that  prevlicted  by  the  G.1-.  SPINNER  program  and  not  within 
75  percent  ot  the  SPINNER  value,  the  SPINNER  value  is  chosen. 

(  omparisons  for  the  Army-Navy  Spinner  test  vehicle  are  made  in  Figure  37.  The  original 
l.MSf  prediction  is  shown  to  depart  consiilerably  from  the  data  and  the  SPINNER  prediction. 


DIMENSIONS  IN  INCHES 


cp  (cal) 


Figure  34.  x<,p  Comparison  for  Air  Slew  Demonstrator 
Forward  of  Body  Midpoint 


Figure  35.  C^i  Comparison  for  Air  Slew  Demonstrator 
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ANGLE  OF  ATTACK  (deg) 


Figure  36.  Modified  Ha.sic  Finner  Aerodynamic  Comparison 
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Damping  Comparison 


riie  revised  LMSC  prediction  improves  the  prediction.  Comparisons  with  data  when  the  center  of 
rotation  is  sliifted  downstream  are  better. 

In  Figure  38,  tlie  L.MSC  model  is  shown  to  predict  the  increased  damping  at  Moo=  1.0  for  the 
basic  tinner  and  follow  the  data  fairly  well  at  the  higher  Mach  numbers.  The  strip  theory  prediction 
tor  lifting  surface  contribution  to  the  pitch  damping  is  seen  to  be  adequate.  As  was  the  case  for 
the  wing  preiliction,  the  strip  theory  prediction  obtained  is  improved  for  higher  aspect  ratios 
and  Mach  numbers.  Unfortunately,  data  fo'-  bodies  with  two  sets  of  lifting  surfaces  were  not 
available. 

Figure  39  shows  that  the  strip  theory  estimate  is  adequate  for  C\p  predictions  at  high  Mach 
numbers.  The  comparison  with  data  and  the  merging  with  the  potential  model  at  lower  supersonic 
Mach  numbers  is  shown  to  be  adequate.  Again  no  data  are  available  for  vehicles  with  two  sets  of 
lifting  surfaces.  As  was  the  case  for  the  wing-alone  static  coefficient  comparisons,  the  accuracy 
obtained  with  potential  theory  is  better  for  higher  aspect  ratio  and  Mach  number. 

TRANSONIC  NOSE  WAVE  DRAG 

Inviscid  pressure  data  are  very  scarce,  particularly  for  the  transonic  case.  Figure  40  shows  the 
comparison  between  limited  data,  the  old  theory,  and  the  current  theory  for  a  M-1  17  bomb  nose, 
rite  numerical  integration  of  pressure  data  was  somewhat  sensitive  since  there  were  not  quite 
enough  pressure  tap  locations.  Flowever,  the  theoretical  prediction,  which  is  based  on  a  full  poten¬ 
tial  solution,  is  much  better  than  the  original  prediction  based  on  the  solution  of  the  full  Euler 
equations.  The  original  Fiuler  computation  did  not  have  enough  grid  resolution  to  accurately 
describe  the  pressure  variation  on  the  noses  of  blunt  pointed  bodies.  The  current  prediction  is  seen 
to  be  an  improvement  over  the  prediction  of  Reference  1. 

Figure  41  compares  the  older  prediction  method,  the  current  method,  and  data  from  Refer¬ 
ence  44.  Again.  ('^|  is  the  forebody  axial  force  coefficient.  The  values  at  M^^  ~  1  labeled  potential 
theory  are  determined  from  supersonic  small  disturbance  theory  for  the  pressure  drag  plus 
Van  Driest's  estimate  of  friction  drag.  Deviations  of  the  supersonic  small  disturbance  values  from 
the  curves  at  M^  =  1 .2  provide  a  measure  of  the  discontinuity  in  the  two  methods  at  M„  =  1.2.  For 
the  blunter  configurations,  a  definite  improvement  is  obtained  by  using  the  new  table  lookup 
procedure.  The  older  method  ncglecteil  the  effect  ot  bluntness. 

TRANSONIC  NORMAL  FORCE  PREDICTION 

As  indicateil  earlier,  the  Nl  AR  envelope  of  body  parameters  is  limited.  The  afterbody  length, 
l.,\.  has  a  maximum  of  .“i  calibers. 

Figure  42  shows  a  comparison  between  the  oliler  computation,  the  Nl-AR  computation,  and 
an  improveil  Nl  AR  computation  with  data.  Fhe  set  ol  bodies  from  Reference  44  with  b-caliber 
afterboilies  is  just  beyond  the  NEAR  envelope  maximum  of  .‘i.  In  general,  the  Cj^j^  prediction  for 
this  set  of  data  is  poor  for  the  NIiAR  computation  anil  improved  by  the  revised  fits  at  the  M^  =  1 .2 
interpolation  point.  Fhe  NFiAR  x^p  prediction,  on  the  other  hand,  is  better  than  the  older  method 
and  considerably  improved  by  the  revised  curve  fit  at  supersonic  Mach  numbers.  As  earlier 
indicated,  the  modified  NEAR  functionali/ation  is  an  input  option  for  less  than  5  calibers. 


Figure  40.  M-1 17  Bomb  Nose  Wave  Drag  Comparison 


Additional  comparisons  for  projectile  shapes  are  given  in  Figures  43  through  45.  All  three 
cases  have  data  with  considerable  scatter.  There  seems  to  be  no  sufficient  justification  for  recom¬ 
mending  any  of  the  methods  over  the  other.  Currently  the  choice  of  algorithm  is  a  program  input 
option.  For  afterbodies  longer  than  5  calibers  the  older  method  is  automatically  used. 

CONCLUDING  RE.MARKS 


The  capabilities  of  an  earlier  aerodynamic  prediction  methodology  have  been  extended  to 
higher  Mach  numbers  and  angles  of  attack.  At  small  angles  of  attack,  the  inviscid  static  aero¬ 
dynamic  prediction  has  been  extended  to  higher  Mach  numbers  (approximately  M^  =  6)  and  the 
pitch  damping  coefficient  is  computed  by  a  new  routine  for  Mach  numbers  from  0  to  approxi¬ 
mately  6  for  body-alone  or  body-tail  configurations.  Use  of  this  routine  is  a  program  option. 
Also  for  small  angles  of  attack,  the  roll  damping  and  pitch  damping  contributions  of  lifting  surfaces 
are  estimated  by  a  strip  theory  at  higher  Mach  numbers  (to  approximately  M^  =6).  This  extends 
the  capability  for  the  roll  damping  computation.  The  capability  for  preiliction  of  the  pitch  damping 
is  extended  to  higher  Mach  numbers  for  conligurations  with  body-canard-tail  configurations  and 
for  body-tail  configuraticjns  when  the  new  pitch  damping  option  is  not  chosen. 

The  transonic  no.se  wave  drag  routine  has  been  replaced.  Significant  improvement  is  obtained 
for  blunted  tangent  ogive  and  blunted  cone  noses.  Also  the  transonic  normal  force  prediction 
routine  has  been  changed.  Use  of  this  routine  is  a  program  option. 

An  empirical  high  angle  of  attack  routine  has  been  adapted  for  body-tail  configurations. 
I  his  routine  is  also  an  input  option.  The  prediction  technique  is  applicable  for  a  limited  envelope 
and  0.8  <  M^.,  <  3  Body-tail  static  aerodynamics  are  predicted  to  45°  angle  of  attack  with 
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Comparison  for  Blunted  Tangent  Ogive-Cylinders 


for  Blunted  Tangent  Ogive-Cylinders  (Continued) 


i 


TRANSONIC  Cm  AND  X^p 


Figure  42.  Cmc,  and  x^p  Comparison  for  Blunted  Tangent  Ogive-Cylinders 
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NORMAL  FORCE  CURVE 
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ij 


Figure  44.  Comparison  of  Theory  and  Test  Data  for  175  mm  XM437  Projectile 
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arbitrary  roll  orientation.  Individual  components  (including  the  body-alone  case)  aerodynamics 
are  predicted  to  180“  for  the  plus  position  roll  orientation. 

A'l  of  the  now  routines  have  been  i..corporated  into  the  existing  code.  Input  logic  and  options 
have  been  improved.  A  detailed  description  of  the  computer  program  and  a  listing  will  be  given  in  a 
later  publication.  Computation  times  depend  on  configuration,  code  option,  and  Mach  number. 
A  computation  time  can  range  from  under  a  second  to  between  30  seconds  and  a  minute  per 
■Mach  number  on  the  CDC  6700. 
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LIST  OF  SYMBOLS 


A 

b 

Ca 

Caf 

^Ap 

^Do 

^DW 

Cg 

Cn, 

^mq  ^  '-  ma 

Cn 

Cnb 

Cnti 

(Cn)tot 

Cp 

Cr 

C, 

D 

D,, 

*BT 

K 

C 

Ca  ’  Lfg  ,  Lp 
M 

Moo 

P 

r 

R 


Lifting  surface  aspect  ratio 
Span  of  two  fins  (excluding  body) 

Distance  along  characteristic  of  the  first  family 

Axial  force  coefficient 

Forebody  axial  force  coefficient 

Pressure  axial  force  coefficient 

Zero  angle-of-attack  drag  coefficient 

Wave  drag 

Roll  moment  coefficient 
Roll  damping  coefficient 
Pitching  moment  coefficient 
Pitch  damping  coefficient 
Normal  force  coefficient 
Body-alone  normal  force  coefficient 

Isolated  ith  fin  (in  presence  of  the  body)  normal  force  coefficient  (cruciform 
configuration) 

Body-tail  total  normal  force  coefficient 

Pressure  coefficient 

Fin  root  chord 

Fin  tip  chord 

Body  reference  diameter 

Mean  body  diameter  near  a  canard  and  tail,  respectively 
Body-tail  normal  force  coefficient  interference 
Hypersonic  similarity  parameter 
Body  length  (cals.) 

Lengths  of  afterbody,  nose,  and  boattail  respectively  (cals.) 

Mach  number  (local) 

Freestream  Mach  number 
Static  pressure 

Cylindrical  cooruinate  (cals.) 

Body  cylindrical  radius  (cals.) 
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Body,  nose  spherical  radius 


Rn 

s  Arc  distance  dong  body 
X  Body  coordinate  (parallel  to  body  axis) 

Xcg  Distance  from  nose  to  moment  reference 
Xj.,  x^^.  Distances  to  fin  apex  for  canard  and  tail,  respectively,  fin  nose 
x^.p  Distance  from  moment  reference  to  center  of  pressure 
Xj^^pij  Body-alone  x^.p  from  nose 
(x^p).j..  Isolated  i"'  fin  (cruciform  configuration)  x^.p  from  nose 
Total  body-tail  x^p  from  nose 
x^^^pi  Interference  x,,,  from  nose 
y  Distance  out  from  a  wing  root 
a  Angle  of  attack 
p  I 

7  C'alorically  perfect  gas  heat  capacity  ratio 
5  Local  body  angle 
Boattail  angle 

X  =  C,/C 

A]  =  Leading  edge  sweep  of  lifting  surface 

M  ==  sin-'  1 

(p  Roll  angle 

Superscripts 

Differentiation  with  respect  to  x 

Subscripts 

a  Asymptotic  value 
w  Wing  or  tail  value 
o  Stagnation  value 

1 .  2  Values  upstream  and  downstream  from  a  body  comer  respectively 
O'  Partial  differential  with  respect  to  a(a  0) 

Reference  area  based  on  D.  Moment  reference  length  based  on  D.  Dynamic  derivative  moment 
coefficients  based  on  (radian  rate)  D 
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